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Abstract 


The  fundamental  fluctuations  of  physical  cross  sections  and  process  rates  in  Quantum 
Electrodynamics  are  known  as  quantum  1/f  effect,  are  limiting  most  high-technology 
applications,  and  are  described  by  the  simple  quantum  1/f  formula  for  the  fractional  spectral 
density,  S5o/a(f)=2aAyfN.  In  spite  of  their  basic  character,  they  are  a  special  case  of  quantum 

1/f  chaos  in  nonlinear  systems,  with  a=e^/hc,  A=(2/37r)(Av/c)2,  Av  the  velocity  change  in 
scattering,  and  N  the  number  of  carriers  used  to  define  the  scattering  cross  section  a.  The 
present  report  provides  for  the  first  time  a  sufficient  criterion  for  the  presence  of  a  1/f 
spectrum  in  an  arbitrary  nonlinear  system,  proves  its  sufficient  character,  and  exemplifies  it 
for  the  current  carriers  in  semiconductors,  for  electrons  in  metals,  for  cars  on  the  highway, 
and  for  the  nonlinear  system  of  particles  and  field  in  electrodynamics.  The  report  gives  first 
principles  analytical  calculation  results  for  the  quantum  1/f  mobility  fluctuations  in  electronic 
materials  such  as  Si  and  GaAs,  derived  on  the  basis  of  the  quantum  1/f  cross  correlation 
formulas  obtained  in  the  previous  AFOSR  Grant  85-0130,  These  quantum  1/f  results  give  the 
conventional  quantum  1/f  noise  as  an  analytical  function  of  temperature,  applied  field,  and  all 
the  physical  parameters  of  the  semiconductor  material,  allowing  for  optimization  of  the 
material  for  any  given  application.  The  quantum  1/f  cross-correlations  derived  by  the  author 
have  been  used  to  recalculate  and  to  graph  the  quantum  1/f  mobility  fluctuations  in  Si  and  GaAs 
samples  as  a  function  of  temperature  and  doping,  in  good  agreement  with  the  measurements  of 
Tacano  in  Japan  and  Hooge  m  the  Netherlands.  The  fractal  dimension  of  band-limited  quantum 
1/f  noise  was  determined  for  the  first  time  and  found  to  coincide  with  the  number  of  octaves 
considered,  in  agreement  with  known  expenmients.  This  indicates  a  quantum  chaos  nature  of  1/f 
noise  in  infrared  detectors,  and  we  expect  this  to  be  valid  in  most  cases  in  semiconductors.  The 
quantum  1/f  theory  was  applied  to  a  quartz  resonator  directly  for  the  first  time,  providing  both 
for  an  explanation  of  the  observed  1/f  frequency  fluctuations  and  for  optimization  means.  The 
practical  application  of  the  quantum  l/f  effect  to  quartz  resonators  in  this  report  allows  us  to 
understand  and  to  extend  the  stability  limits  of  quartz  resonators.  A  fundamental  breakthrough 
was  performed  through  the  first  direct  derivation  of  the  coherent  quantum  1/f  effect  from  a 
special  quantum-electrodynamic  propagator,  and  from  the  author's  general  sufficient  1/f  chaos 
criterion  presented  in  this  report.  Finally,  a  suggestion  is  presented,  which  provides  an 
explanation  of  the  discrepancies  which  have  been  noticed  in  the  coiiector  current  1/f  noise  of 
bipolar  junction  transistors.  Our  new  formula  for  collector  1/f  noise  in  ultrasmall  BJT's  was 
found  to  agree  reasonably  with  the  experiment. 
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I.  INTRODUCTION 


Progress  has  been  achieved  through  this  grant  in  the  study  of  nonlinear  systems  which 
generate  chaotic  1/f  fluctuations,  in  the  application  of  the  Quantum  1/f  Theory  [1-13]  to 
various  materials  used  in  small  and  ultrasmall  electronic  devices,  in  the  application  of  the 
Quantum  1/f  Theory  to  electronic  devices  and  in  the  fundamental  Quantum  1/f  Theory  itself. 

The  close  similarity  of  the  classical  and  quantum  1/f  theories  and  the  initial 
development  of  the  quantum  1/f  theory  by  the  author  out  of  his  efforts  to  quantize  his  classical 
turbulence  theory,  have  led  to  sustained  efforts  of  the  author  aimed  to  integrate  all  his  theories 
as  various  forms  or  realizations  of  a  fundamental  notion  of  chaos  in  nonlinear  systems.  During 
this  grant  period,  these  efforts  finally  beared  fruits.  A  general  sufficient  criterion  was 
formulated,  allowing  to  identify  the  nonlinear  systems  which  exhibit  1/f  spectra.  This 
criterion  is  presented  in  Sec.  II  below.  It  displays  for  the  first  time  the  general  mathematical 
principle  causing  the  ubiquitous  1/f  spectrum  in  nature,  science  and  technology.  It  is  followed 
in  Sec.  Ill  by  examples,  in  which  the  new  criterion  is  applied  to  the  classical  and  quantum 
mechanical  forms  of  the  author's  1/f  noise  theory.  These  five  examples  clarify  the  physical 
meaning  of  the  new  criterion. 

Infrared  detectors  have  provided  a  testing  ground  for  the  quantum  1/f  theory  in  the 
eighties,  and  have  been  the  subject  of  several  quantum  1/f  studies.  In  Sec.  IV  we  present  a 
practical  application  of  the  quantum  1/f  theory  to  infrared  detectors  in  the  presence  of 
radiation,  in  particular,  a  proof  of  the  absence  of  quantum  1/f  noise  in  the  process  of  carrier 
photogeneration  in  semiconductors.  Quantum  1/f  noise  is  nevertheless  predicted  by  the 
quantum  1/f  theory  and  is  present  both  in  junction  and  MIS  photodetectors,  in  the  dark  current 
and  in  the  resistance  of  the  homogeneous  part  of  the  semiconductors  which  form  the  detector. 

Quantum  1/f  noise  was  introduced  as  the  quantum  form  of  the  classical  turbulence 
process  which  is  a  form  of  chaos.  Quantum  1/f  noise  is  therefore  a  fundamental  and  universal 
form  of  quantum  chaos.  The  first  calculation  of  the  fractal  dimension  of  the  quantum  1/f  noise 
process  presented  in  Sec.  V. 

The  quantum  1/f  theory  has  inspired  the  explanation  of  the  origin  of  the  famous  Q''^  law 
of  flicker  of  frequency  in  quartz  crystal  resonators  by  the  author,  as  early  as  1978.  He  has 
now  developed  a  new  method  allowing  for  direct  application  of  the  quantum  1/f  principle  to 
quartz  resonators  and  to  other  piezoelectric  and  ferroelectric  systems,  which  is  presented  in 
Sec.  VI.  That  section  includes  the  successful  practical  application  of  the  quantum  1/f  theory  to 
quartz  resonators  and  the  verification  of  the  results  at  NIST-Boulder  with  the  cooperation  of 
Fred  Walls.  The  quantum  1/f  effect  appears  in  the  process  rate  of  phonon  removal  from  the 
main  oscillator  mode  of  the  quartz  crystal,  because  there  is  bremsstrahlung  associated  with  this 
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process  in  any  piezoelectric  crystal.  Indeed,  as  we  shall  see  in  Sec.  VI,  the  process  will 
suddenly  change  not  only  the  velocity  of  the  lattice  vibration,  but  also  the  rate  of  chance  of  the 
polarization  of  the  crystal,  which  in  turn  is  equivalent  to  sudden  electric  current  change. 

We  have  tried  to  improve  the  application  of  quantum  1/f  theory  to  the  collector  noise  of 
bipolar  transistors.  This  short  calculation  is  presented  in  Sec.  VII  which  also  provides  a 
comparison  of  our  improvements  in  the  quantum  1/f  collector  noise  formula  in  ultrasmall  BJTs 
with  the  experimental  evidence. 

For  the  practical  application  of  the  Quantum  1/f  Theory  it  is  necessary  to  derive  the 
quantum  1/f  fluctuations  of  various  kinetic  (transport)  coefficients  which  characterize  the 
materials  used  in  electronic  and  microelectronic  applications  (mobility,  diffusion  coefficient), 
from  the  author's  fundamental  quantum  1/f  formula.  The  latter  is  applicable  only  to  cross 
sections  and  rates  of  elementary  processes.  Most  important  is  the  calculation  of  mobility 
fluctuations  in  Si,  GaAs  and  Hgi-xCdxTe.  An  earlier  calculation  (Kousik,  Van  Vliet,  Handel, 
1985)  of  mobility  fluctuations  in  Si  and  GaAs  is  replaced  in  Sec.  VIII  by  a  more  rigorous 
calculation,  based  on  the  new  quantum  1/f  cross-correlations,  developed  under  the  previous 
AFOSR  Grant,  and  presented  in  the  Final  Technical  Report  AFOSR  -85-0130.  The  new 
calculation  yields  increased  1/f  noise,  and  is  in  vey  good  agreement  with  the  experiment.  The 
results  give  the  quantum  1/f  noise  present  in  various  carrier  scattering  mechanisms  as  a 
function  of  temperature  and  doping  level. 

The  progress  achieved  in  the  fundamental  quantum  1/f  theory  includes  the  direct 
derivation  of  the  Coherent  Quantum  1/f  Effect  from  a  special  quantum-electrodynamical 
propagator  known  since  1975,  and  from  the  sufficient  1/f  criterion  introduced  in  the  1992 
(Third)  Annual  Report  (Sec.  II).  The  physical  meaning  of  the  coherent  quantum  1/f  effect  is 
displayed  by  this  new  independent  derivation  as  a  consequence  of  the  indefinite  energy  of  a 
physical  electron. 

At  the  same  time,  many  new  contributions  to  quantum  1/f  theory  and  experiment  were 
published  by  other  workers  in  the  field,  considerably  advancing  the  field  of  infra-quantum 
physics  and  quantum  1/f  noise  in  high-technology  applications.  These  new  contributions,  as 
well  as  new  PhD  thesis  work  in  this  field  and  contributions  presented  at  the  “4^^  Conference  on 
Quantum  1/f  Noise  and  other  Low-Frequency  Fluctuations"  (Minneapolis  1990),  at  the  5^^  (St. 
Louis,  1992)  and  at  the  X,  XI  and  XII  Int.  Conf.  on  Noise  in  Physical  Systems  (Budapest,  1989, 
Kyoto,  1991,  St.  Louis,  1993)  are  all  included  in  the  updated  General  Quantum  1/f 
Bibliography  appended  to  this  Report.  My  collaborators  have  been  A.L.  Forst-Chung,  E. 

Bernard!,  T.  Chung,  L.M.N.  Sastri,  X.  Hu,  Jian  Xu,  Y.  Zhang,  and  M.  Leong. 
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II.  SUFFICIENT  CRITERION  FOR  1/F  NOISE  IN  CHAOTIC  NONUNEAR  SYSTEMS 


Consider  a  n-dimensional  nonlinear  system  described  in  terms  of  the  dimensionless  function 
Y(x,t)  by  the  m^^  order  nonlinear  dynamical  equation 

dY/dt  +  F(x,  Y,  dY/dx-j  ...dY/dXn,  d2Y/dxi2 . dn^Y/dXn"^)  =  0  (  1  ) 

If 

F[Xx,  Y,  dY/(>idx-|)...dY/(?LdXn).  d2Y/{Xdxi)2 . d^^ Y/(XdXn)^] 

=  ?iPF(x,  Y,  dY/dxi  ...dY/dxn,  d2Y/dxi  2 . d^ Y/dXn'^)  (  2  ) 

for  any  real  number  X,  Eq.  (1)  is  said  to  be  homogeneous.  Performing  a  Fourier  transformation 
with  respect  to  the  vector  x,  we  get  in  terms  of  the  Fourier-transformed  wavevector  k  the 
nonlinear  integro-differential  equation 

dy(k,t)/dt  G[k,  y(k,t),  ki  y(k,t)...kny(k,t).  k-|2y(k,t) . k^  y  (k.t)]  =  0,  (3) 

where  y{k,t)  is  the  Fourier  transform  of  Y(x,t).  Due  to  Eq.  (2),  G  satisfies  the  relation 

G[Xk,  y,  Xkiy...Xkny.  (Xki)2y . (Xkn)^y] 

=  XPG[k,  y,  kiy...kny,  ki2y . kn^^y].  (  4  ) 


Eq.  (3)  can  thus  be  rewritten  in  the  form 

dy/d(t//.P)  +  G[Xk,  y,  Xk^y...XknY,  (Xki)2y . (Xkn)f^y]  =  0,  (5) 

Taking  >.=1/k,  where  k=|k|=(x-i 2+ __.+Xn2)l/2_  sncj  setting  kPt=z,  we  notice  that  k  has  been 
eliminated  from  the  dynamical  equation,  and  only  k/k  is  left.  This  means  that  there  is  no 
privileged  scale  left  for  the  system  in  x  or  k  space,  other  than  the  scale  defined  by  the  given 
time  t,  and  expressed  by  the  dependence  on  z.  We  call  this  property  of  the  dynamical  system 
"sliding-scale  invariance". 

In  certain  conditions,  instabilities  of  a  solution  of  Eq.  (1)  may  generate  chaos,  or  turbulence. 
In  a  sufficiently  large  system  described  by  the  local  dynamical  equation  (1),  in  which  the 
boundary  conditions  become  immaterial,  homogeneous,  isotropic  turbulence,  (chaos)  can  be 
obtained,  with  a  spectral  density  determined  only  by  Eq.  (1).  The  stationary  autocorrelation 
function  A(x)  is  defined  as  an  average  over  the  turbulent  ensemble 
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A(x)  =  <Y(x,l)Y(x,t+x)>  =  |<y(k,t)y(k,t+T)>d'^k  =  Ju(k,z)d^k 


(6) 


Here  we  have  introduced  the  scalar 

u(k,z)  =  <y(k,t)y(k,t+T)  >  (7) 

of  homogeneous,  isotropic  chaos  (turbulence),  which  depends  only  on  |k|  and  z=kPx.  All 
integrals  are  from  minus  infinity  to  plus  infinity.  The  chain  of  integro-differential  equations 
for  the  correlation  functions  of  any  order  obeys  the  same  sliding-scale  invariance  which  we 
have  noticed  in  the  fundamental  dynamical  equation  above.  Therefore,  in  isotropic, 
homogeneous,  conditions,  u  can  only  depend  on  k  and  z.  Furthermore,  the  direct  dependence  on  k 
must  reflect  this  sliding-scale  invariance,  and  is  therefore  of  the  form 

u(k,z)  =  k'^v(z).  (  8  ) 

Indeed,  only  this  form  insures  that  u(k,z)d^k  and  therefore  also  the  corresponding  integrals  and 
multiple  convolutions  in  k  space  have  the  necessary  sliding-scale  invariance. 

According  to  the  Wiener-Khintchine  theorem,  the  spectral  density  is  the  Fourier-transform 
of  A(t), 

Sy(f)  =  je2^iftA(x)dx  =  (l/f)  Je2Ji‘t'|k-nv(z)dnk’dt'  =  C/f,  (9) 

where  we  have  set  fx=t',  k^=fk''^,  z=k^x=k'^t',  and  the  integral 

C  =  Ie2rrit'|k-nv(z)dnk'dt'  =  le2j^ifjk-nv(k”n)dnk"dt'  (10) 

is  independent  of  f.  We  have  defined  the  vector  k"=t’'''''^  k. 

In  conclusion,  we  have  shown  that  if  the  equation 

dY/dt  +  F(x,  Y,  dY/dxi...dY/dxn.  . df^Y/dxnf^)  =  0 

with 

F[Xx,  Y,  dY/(Xdxi)...dY/(;idx)p,  d2Y/(Xdxi)2 . d^ Y/(5idXn)^]  =  ?iPF(x,  Y, 

dY/dx-i  ...dY/dxp,  d2Y/dxi2 . dn^Y/dxp^^) 

admits,  in  the  limit  of  weak  dissipation,  quasistationary  homogeneous  isotropic  chaotic 
(turbulent)  solutions  which  are  practically  independent  of  the  nature  of  the  instabilities  or 


(11  ) 

(12) 
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bifurcations  (or  even  stirring  forces)  which  have  caused  the  chaotic  state,  the  corresponding 
spectral  density  must  be  proportional  to  1/f. 

We  note  that  the  solution  (7)  leads  to  a  weak  (logarithmic)  divergence  of  the  integral  over  k 
in  the  last  form  of  Eq.  (6)  and  in  Eq.  (9).  This  seems  to  contradict  the  fact  that  in  practical 
applications  the  autocorrelation  function  A(t)  is  finite,  and  its  value  at  t=0  is  usually  given  in 
the  problem  at  hand.  However,  in  practice  one  never  deals  with  an  infinite  volume,  and  the 
physical  wave-vectors  are  also  limited.  For  instance  in  fluid  dynamics,  wave  vectors  exceeding 
the  reciprocal  average  distance  between  neighboring  fluid  molecules  correspond  to  thermal 
motions,  and  are  therefore  no  longer  meaningful  for  the  hydrodynamic  treatment.  Due  to  its 
logarithmic  character,  the  divergnce  is  thus  without  practical  importance.  Nevertheless,  for  a 
given  level  of  chaos  A(0),  we  can  construct  an  approximate  solution 

u(k,z)=kr-nv(z),  (13) 

with  0<r«1,  which  avoids  the  divergence  at  k=0.  To  get  the  correct  chaos  level  with  k<kQ,  ko 
being  an  upper  cutoff,  we  set 

u  =  r[A(0)/v(0)]ko''’k'''^v(2) .  (14) 

This  yields  for  t=0  the  result  A(0)  when  we  integrate  over  d^k  with  an  upper  limit  kQ.  We 
notice  that  r  is  present  both  as  a  general  factor,  and  as  a  small  defect  in  the  exponent  of  k.  This 
is  a  general  feature,  present  both  in  classical  and  quantum  nonlinear  systems  with  1/f  noise.  In 
the  limit  r->0,  the  approximate  solution  tends  to  become  exact.  In  the  classical  homogeneous, 
isotropic,  turbulence  theory,  r  can  be  arbitrarily  small,  while  in  the  quantum  1/f  theory 
(quantum  electrodynamics),  r  -  aA  =  (2a/3r)(Av/c)2  «  1/137  is  the  well-defined  infrared 
exponent  of  the  process,  with  ci-e2/hc=in37,  as  we  know  from  the  theory  of  infrared 
radiative  corrections.  This  allows  us  to  formulate  the  justification  of  our  criterion  in  a  way 
used  by  the  author  for  the  first  time  in  19S1  in  an  unpublished  paper  which  was  reviewed  in 
that  year  by  T.  Musha  and  C.M.  Van  Vliet;  Homogeneous  nonlinear  systems  require  the  1/f 
spectrum  because  this  is  the  only  one  which  is  idempotent  with  respect  to  the  self-convolution 
operation  in  the  limit  of  very  small  r  values.  Indeed,  the  result  of  the  multiple  self¬ 
convolutions  introduced  by  the  Fourier  transformation  of  a  nonlinear  term  is  again  1/fi'''  with 
a  very  small  r,  if  the  input  was  a  1/f^'’’  with  a  very  small  r. 
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III.  EXAMPLES  ILLUSTRATING  THE  APPLICATION  OF  THE  GENERAL  SUFFICIENT 
CRITERION  FOR  1/f  NOISE  IN  CHAOTIC  NONUNEAR  SYSTEMS 


In  spite  of  the  practical  success  of  our  quantum  1/f  theory  in  explaining  electronic  1/f 
noise  in  most  high-tech  devices,  and  in  spite  of  the  conceptual  success  of  our  earlier  classical 
turbulence  approach  to  1/f  noise,  the  question  about  the  ultimate  origin  of  nature's  omnipresent 
1/f  spectra  remained  unanswered.  During  the  last  three  decades,  we  have  claimed  repeatedly 
that  nonlinearity  is  a  general  cause  of  1/f  noise.  Our  new  criterion  proves  that  nonlinearity 
always  leads  to  a  1/f  spectrum  if  homogeneity  is  also  present  in  the  equation(s)  of  motion.  We 
can  present  this  criterion,  derived  in  the  preceeding  section,  in  a  more  general  form. 
Specifically,  let  the  system  be  described  in  terms  of  the  dimensionless  vector  function  Y(x,t) 
by  the  m'*^  order  nonlinear  system  of  differential  equation 

cr>[t,  X,  Y,  aY/3t,  3Y/axi...aY/3xn,  a2Y/at2,  a2Y/axi2 . a^^Y/axp^^]  =  o  ( 1  ■) 

where  the  vector  function  <I)  may  be  nonlinear  in  any  of  its  arguments.  If  a  number  0  exists 
such  that  Eq.  (1')  implies 

ci)[xet,  J.X,  Y,  aY/xeat,  aY/xaxi  ...aY/xaxp,  a2Y/x26at2, 

a2Y/x2axi2 . af^Y/xniaxn^)  =  0  (2’) 

for  any  real  number  X,  the  power  spectral  density  of  any  chaotic  solution  for  the  vector  function 
Y  defined  by  Eq.  (1')  is  proportional  to  1/f. 

In  conclusion,  nonlinearity  +  homogeneity  =  1/f  noise.  The  ultimate  cause  of  the 
ubiquitous  1/f  noise  in  nature  is  the  omnipresence  of  nonlinearities  (no  matter  how  weak)  and 
homogeneity.  The  latter  is  finally  related  to  rotational  invariance  and  to  the  isotropy  of  space. 
All  our  four  specific  theories  of  1/f  chaos  in  nonlinear  systems  are  just  special  cases  to  which 
our  criterion  is  applicable.  They  include  our  magneto-plasma  theory  of  turbulence  in  intrinsic 
symmetric  semiconductors  (1966.  (14-16)).  our  similar  theory  for  metals  (1971,  [17]), 
the  quantum  1/f  theory  (pure  quantum  electrodynamics,  1975,  [1-13]),  and  the  theory  of 
Musha's  traffic  turbulence  (1989).  A  fifth,  negative,  application  example  concerns  a  one¬ 
dimensional  crystal,  i.e.,  a  chain  of  atoms  with  slightly  anharmonic  interaction  potentials,  in 
which  the  criterion  indicates  the  absence  of  a  lattice-dynamical  1/f  spectrum,  in  agreement 
with  computer  simulations  by  Musha  and  with  analytical  calculations  by  Koch  and  Wolf. 

The  general  criterion  developed  in  the  preceeding  section  will  now  be  illustrated  on  the 
basis  of  the  examples  mentioned  above. 
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111.1  Classical  Turbulence  Theory  for  the  Current 
Carriers  in  Semiconductors 


In  the  case  of  homogeneous,  isotropic  turbulence  [14-17]  caused  in  the  electron-hole 
plasma  of  an  infinite  sample  of  a  symmetric  intrinsic  semiconductor  by  dynamical  instabilities 


of  any  kind,  we  start  from  the  equations 

VV+  =  (e/2c)v-xB  -  (1/n)VP,  (15) 

vv  =  2e[E  v  +  xB/c]  -  (2/n)V(Pp  -  Pn),  (16) 

V-v+=  0  (n  =  const),  (17) 

V  xE  =  -(1/c)aB/at,  (18) 

V  xB  =  2Tcenv/c,  (19) 

VB  =  0.  (20) 


Here  n  is  the  total  carrier  concentration  including  an  equal  number  of  electrons  and 
holes,  v/e  their  reciprocal  mobility  assumed  to  be  the  same  for  electrons  and  holes,  Pn  and  Pp 
the  partial  pressures  of  electrons  and  holes,  P  the  total  carrier  pressure,  2V''-  v'  the  sum  and 
the  difference  of  the  carrier  drift  velocities.  Inertial  terms  proportional  to  the  effective 
masses  of  the  carriers,  as  well  as  electrostatic  terms  and  compressibility  terms  have  been 
neglected  here  in  a  consistent  way  [14-17],  because  we  are  interested  in  the  low-frequency 
domain  only.  Although  we  do  not  work  this  out  here,  this  system  of  equations  can  be  shown  to 


admit  an  energy  theorem.  Performing  a  Fourier  expansion,  we  obtain 

vv'''(k)  =  (e/2c)Z  k’v'(k')xB  (k-k')  -  (i/n)kP(k),  (21) 

vv-(k)  =2e[E(k)-t-  Z  k'V'-(k')xB  (k-k')/c]-(2i/n)k(Pp-Pn),  (22) 

k-v  +  (k)  =  0.  (23) 

ik  xE(k)  =  -(1/c)aB(k)/at.  (24) 


ik  xB(k)  =  2Tce n v(k)/c, 


(25) 


k-B(k)  =  0,  (26) 

Substituting  E  from  Eq.  (22)  into  Eq.  (24),  we  obtain  with  Eq.  (25) 

3B(k)/at  +  4k2B(k)  =  ikx  X  k*  v  +  (k')  xB(k-k’),  (27) 


where  ^  =  c^vMTtne^.  Eqs.  (21)  and  (25)  yield 

v+(k)  =  (i/4;tvn)  Z  k"{  B  (k")[k"B  (k-k”)]  -  k"[B  (k")B  (k-k")] 
-(k/k2)(i  .5k,o)[k"B  (k-k")  kB(k")  -  B  (k")B  (k-k")  (k-k")]}  . 


(28) 


Substituting  this  into  Eq.  (27),  we  obtain  the  fundamental  dynamical  field-equation  of 
turbulence  in  the  electron-hole  plasma  of  a  symmetric  intrinsic  semiconductor 


abp(k,t)/at  +  pk2bp(k,t)  =  Z  k'k"  bj(k-k’,t)  b|(k",t)  bm  (k'-k",t) 
•(kj5ps-ks5p])[ks5|m-km"5ls  +  (ksVk'2)(i  .5k',o)(km  "ki’-k'-k-Sim)]. 


(29) 


in  terms  of  b  =  B/V2:tvn.  This  dynamical  equation  has  the  form  of  Eq.  (3),  with  p=2  in  Eq. 

(4),  and  with  G  defined  as  the  r.h.s.  minus  the  term  in  k2  on  the  l.h.s..  Our  sufficient  criterion 
thus  tells  us  that  this  nonlinear  system  will  yield  a  1/f  spectrum.  We  present  below  the  direct 
derivation  for  this  example. 

Multiplying  Eq.  (29)  with  ba*(k,t-T)  and  taking  the  average  over  a  statistical 
ensemble  which  represents  our  notion  of  stationary  turbulence,  we  obtain  in  quasi-stationary 
conditions 


5wap(k,x)/aT  -I-  pk2wa3(k,-c) 

=  Zk-k"<ba*(k,t-x)bj(k-k’,t)b|(k",t)bm(k'-k",t)>Rjmlp, 


1  1 


(30) 


with  Wap(k,T)  =  (L/27:)3<ba*(k,t-':)bp(k,t)>,  where  L  is  the  edge  of  the  cubic  normalization 
box,  and 

Rjmip  =(kj5ps-ks5pj)[ks5|m-km"5|s  +(ksVk-2)(i -5k',o)(km  ”kr-k'-k"5|m)].  (31  ) 

Multiplying  Eq.  (29)  with  more  magnetic  field  components  and  averaging,  we  obtain 
equations  connecting  the  fourth-order  correlation  tensor  to  the  sixth-order  tensor,  and  so  on 
[14-17].  To  end  this  infinite  chain  of  equations  for  the  correlation  functions,  we  make  a 
quasinormality  assumption  which  expresses  the  fourth-order  moment  appearing  in  Eq.  (30) 
according  to  the  scheme 

<ABCD>  =  <AB>  <CD>  +  <AC>  <BD>  +  <AD>  <BC>,  (32) 


valid  if  the  four  field  components  would  approximate  a  joint  normal  distribution.  This 
approximation  does  not  alter  the  homogeneity  of  the  system,  which  ultimately  causes  the  1/f 
spectrum.  This  approximation  yields  the  closed  equation 

3wap(k,x)/aT  +  iik2wap(k,x) 

=  (2;r/L)3wa](k.x)Ik'W|m(k’.0)Rjmlp(k,k-).  (33) 


Isotropic  turbulence  requires  w^p  =  Ai(k)5ap  +  A2(k)kockp,  with  coefficients  Ai  and  A2 
related  through  Eq.  (26),  yielding 

Wap(k,x)=  (1/2)[5ap  -  kakp/k2]u(k.x).  (34) 


where  u(k,x)  =  SaWaa(k,x).  Therefore,  the  scalar  correlation  function  u(k,x)  satisfies  the 
dynamical  equation  of  homogeneous,  isotropic,  stationary  turbulence 


3v(k,x) 

a|x|  ^ 


v(k,x)  =  -  ^(k.x 


.  fd^k-  k2  +  k  k’r .  /k  k\oi  ..  . 


(35) 


where  v(k,x)Hk'3u(k,x).  and  x  =  nxk2  js  a  dimensionless  variable  replacing  x.  We  convince 
ourselves  that  the  integral  is  independent  of  k.  provided  v(k,x)  does  not  depend  on  its  first 
argument,  by  setting  k'/k^K.  This  yields  a  solution.  However,  with  v=e’^l^l  we  get  a 
logarithmic  divergence  at  k=0.  We  look  for  an  exact  solution  of  the  form  [14-17] 


v(k,x)  =  or  u(k,x)  =  (h/k^-^) 


(36) 


where  m(k)  is  very  close  to  a  constant,  almost  independent  of  k,  and  h  is  a  constant  proportional 
to  the  intensity  of  the  turbulence,  or  the  turbulence  level.  Substituting  this  into  Eq.  (35)  and 
performing  the  integration,  we  obtain  a  finite  result  only  for  0<e<2: 

m(k)  =  1  +  h  r(e)k^,  with  r(e)  =  [27t2cotan(en/2)]/[(1 -e2)(3-e2)].  (37) 

We  notice  that  m(k)  is  indeed  practically  constant  when  0<e«1  is  very  small,  arbitrarily 
small.  The  value  e=0  leads  to  a  logarithmic  divergence,  but  we  can  set  e=0  for  practical 
purposes. 

The  spectral  density  corresponding  to  Eq.  (36)  with  £=0  is 


Wap(“)  =  (1/rt)  Jcoscoidt  Jwo,p(k,x)d2k 


oc 

-  jSaU  Jk^dkldi  u(k,i)C0SMt  . 


(38) 


This  is  a  1/f  spectrum.  At  the  low  frequency  end  we  do  not  get  a  divergent  spectral  integral, 
because  the  more  exact  form  of  the  spectrum  with  a  finite  small  e«1  is  [14-17] 

oo  oo 

r mk'*'*'^dk  _ 1 _ 

Jco^-t-m^k'^  u  ^ 2  J  i+x^  ’  (29) 

0  0 

which  is  proportional  to  f^^2-1  n  jg  interesting  to  note  that  for  e«1  cotan  E7r/2  =  2/e7:  in  Eq. 
(37),  and  that  the  value  of  e  calculated  from  Eq.  (37)  is  therefore  proportional  to  h,  or  to  the 
intensity  of  the  turbulence.  This  feature  of  the  classical  theory  [14-17]  is  expressed  with 
fascinating  clarity  in  the  quantum  form  of  the  theory,  where  e  is  replaced  by  2aA  which  also 
appears  as  a  intensity  factor  multiplying  the  quantum  1/f  noise. 

The  essential  element  which  led  to  the  1/f  spectrum  in  the  classical  turbulence  theory  is 
the  nonlinearity  of  the  equations  of  motion,  caused  by  the  reaction  of  the  electric  currents  back 
on  themselves  via  the  generated  electromagnetic  field.  The  same  feedback  reaction,  via  the 
electromagnetic  field,  also  caused  the  nonlinearity  in  the  quantum  1/f  theory,  and  in  QED  in 
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general,  leading  in  the  same  way  to  an  identical  1/f  spectrum,  this  time  with  a  physically  more 
meaningful  e=2aA.  This  nonlinearity  induces  the  coupling  between  various  scales  of  turbulence 
and  leads  to  the  dynamical  equilibrium  between  eddies  of  all  sizes,  expressed  by  the  1/f 
spectrum,  in  the  e=0,  or  aA=0,  limit,  this  dynamical  equilibrium  assumes  both  for  the 
quantum  case  and  for  homogeneous,  isotropic,  turbulence  in  the  unbounded  semiconductor 
sample  the  simplest  form,  characterized  by  scale-homogeneity,  or  scale  invariance.  Indeed, 
replacing  for  e=0  in  Eqs.  (35)  and  (36)  k  and  k'  by  Xk  and  Xk',  while  also  replacing  x  by  x/X^, 
(or  0)  by  Eq.  (35)  is  not  affected,  and  X  drops  out.  We  conclude  that  in  the  weak 
turbulence  limit  (e=0)  we  obtain  perfect  self-similaritv  of  the  turbulence  proceee  at  all  scales 
in  space  and  time,  classically  and  quantum-mechanically.  The  implied  scale  invariance  is 
caused  by  the  absence  of  any  characteristic  length  or  time  scale,  or  by  the  presence  of  a  sliding 
scale.  Indeed,  the  frequency  scale  pk^  is  a  function  of  the  size  of  the  eddies,  given  by  the  wave 
number  k  which  can  have  any  value.  The  actual  frequency  and  wave-number  spectra  are  closely 
related  fractals,  but  in  the  weak-turbulence  limit  they  approach  an  exact  1/f  and  l/k^ 
spectrum  respectively.  In  fact,  we  are  here  understanding  the  nonlinear  dynamics  which 
shapes  this  fractal  for  the  first  time. 

III. 2  Turbulence  Theory  for  Drude  Electrons  in  Metals 

Our  classical  turbulence  theory  can  be  extended  to  the  case  of  metals  or  degenerate 
extrinsic  semiconductors  [17]  in  the  Drude  free  electron  gas  model.  The  system  of  integro- 
differential  equations  is  quite  different, 


vv'*'  =  -eE  -  (e/c)vxB  -  (1/n)VP 

(40) 

VxE  =  -(1/c)aB/af 

(41  ) 

VxB  =  -(47ren/c)v 

(42) 

V-B  =  0, 

(43) 

and  leads  to  a  third-order  nonlinearity  [17]  in  the  resulting  closed  equation  of  motion,  or 
nonlinear  field-equation,  which  replaces  Eq.  (28): 

aB(k,t)/at  +  vk2B(k.t)  =  -(c/47ine)kxJd3k'B(k-k‘,t)x[k'xB(k',t)]. 
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(44) 


This  is  again  in  the  form  of  Eq.  (3),  with  p=2  in  Eq.  (4).  We  thus  expecta  1/f  spectrum  in  this 
system  as  well.  This  time  we  only  sketch  the  derivation.  The  corresponding  infinite  chain  of 
equations  for  the  correlation  tensors  now  goes  in  steps  of  one.  As  was  shown  above,  it  went  in 
steps  of  two  for  semiconductors.  The  third-order  correlation  can  be  eliminated  between  the 
first  and  second  equations  in  the  chain.  The  resulting  dynamical  equation  [17]  for  homogeneous, 
isotropic,  stationary  turbulence,  which  replaces  Eq.  (35),  with  the  same  notations,  using  03  as 
the  unit  vector  of  the  third  axis,  is 


3^v(k,x) 

3x2 


-(-  v(k,x) 


=  -  ^2  ]v(Kk,K2x)v(k|e3  +  K|,xle3  +  K|2). 


ea+K 


(45) 


This  also  admits,  in  the  e=0  limit  of  weak  turbulence,  a  solution  v(k,x)  which  does  not 
depend  on  the  first  argument,  and  u(k,x)  =  k‘3e'^'^(^),  this  time  with  an  x-dependent  m. 

With  the  change  of  variables  t=a)T  and  k'=k/Vco  in  the  second  (middle,  involving  u)  form  of  Eq. 
(38),  X  remains  invariant,  and  a  factor  1/co  will  appear  in  front  of  the  integrals  which 
themselves  will  just  yield  a  constant  factor  independent  of  co,  t  or  k.  We  thus  obtain  again  a 
universal  1/f  spectrum.  As  is  shown  in  detail  elsewhere  [17],  this  1/f  spectrum  is  expressed 
in  the  corresponding  current  and  voltage  fluctuations  which  can  be  observed  in  the  semiconduc¬ 
tor  or  metallic  medium.  We  conclude  that  the  1/f  spectrum  is  a  general  property  of  electrically 
conducting  systems  in  interaction  with  the  electromagnetic  field,  a  property  which  is  caused  by 
the  nonlinearity  of  the  system  of  carriers  and  field  in  mutual  interaction  due  to  the  absence  of  a 
characteristic  scale  in  the  nonlinear  equation  of  motion,  and  which  finds  its  clearest  expression 
in  the  Quantum  1/f  Effect. 


111.3  Theory  of  Highway  Traffic  Fluctuations 

Musha  and  Higuchi  [18,19]  discovered  the  1/f  spectrum  of  highway  traffic  fluctuations 
empirically  in  1977,  also  developing  a  model  based  on  a  postulated  linear  dependence  of  the 
average  traffic  speed  v  on  the  linear  concentration  of  cars  n(x)  on  the  road,  v=vo(1-n/ns). 

The  model  was  recognized  to  be  similar  to  Burger's  model  of  turbulence,  and  was  simulated 
numerically  leading  to  a  1/f-like  spectrum  at  low  frequencies. 

The  present  paper  develops  a  statistical  turbulence  theory  for  the  Musha  model,  showing 
how  in  the  low  wave  number  and  low  frequency  limit  a  universal  1/f  spectrum  emerges.  In  this 
limit,  this  is  a  particular  case  of  the  author's  sliding-scale-invariant  class  of  nonlinear 


systems,  all  characterized  by  the  presence  of  a  universal  1/f  spectrum.  The  author's  classical 
and  quantum  1/f  theory  is  another  example  in  the  same  class. 

Musha  writes  the  traffic  current  J  as  a  sum  of  drift  and  diffusion  currents 

J  =  nv  -  D3n/3x  =  nvo  -  (vo/ns)n2  -  D3n/3x,  (4  6) 

where  D  is  a  "diffusion  coefficient".  The  equation  of  continuity  is  [8,9] 

0  =  3n/3t  +  3J/3X  =  (3/3t  +  Vo3/3x)n  -  2(vo/ns)n3n/3x  -D32n/3x2.  (47) 

This  is  Musha's  fundamental  equation  of  traffic  dynamics,  which  also  was  written  [18]  in  a 
system  of  reference  defined  by 

X'  =  -X  +  vot,  t'  =  t  (4  8) 

in  the  final  form  [18] 

3n/3t'  +  an3n/3x'  =  D32n/3x'2,  (4  9) 

where  a  s  2vo/ns. 

We  normalize  the  concentration  to  ns/2,  the  coordinate  along  the  road  x  to  D/vq,  and  the 
time  to  D/vq^,  thereby  obtaining  the  dimensionless  form  of  the  fundamental  traffic-dynamical 
equation 

3n(x,t)/3t  +  n(x,t)3n(x.t)/3x  >  32n(x.t)/3x2.  (50) 

where  we  did  not  bother  to  change  the  notation,  and  returned  to  the  original  unprimed  notation. 
Expanding  in  a  Fourier  series  over  the  interval  L.  we  obtain 

3n/3t  +  iZk’k'n(k')n(k-k‘)  =  -k^nfk.t).  (51) 

Comparing  Eq.  (50)  with  Eqs.  (1),(2),  or  Eq.  (51)  with  Eqs.  (3), (4),  we  notice  that 
this  time  our  criterion  is  not  satisfied,  due  to  the  r.h.s.  term.  However,  in  the  limit  of  small  k 
that  term  becomes  negligible,  and  thus,  our  criterion  becomes  applicable,  and  we  should  get  a 
1/f  spectrum  in  the  low  frequency  limit.  This  would  correspond  to  p=1  in  Eqs.  (2)  and  (3). 
We  derive  this  1/f  spectrum  below. 


Defining  the  autocorrelation  function 


A(^,t)  =<n*(x,t)n(x+^,t+x)>  =Zi<<n*(k,t)n{k,t+x)>e‘*^^  =ZkU(k,x)e'*^^,  (52) 

as  a  turbulent  ensemble  average  <...>,  in  quasistationary  homogeneous  conditions  we  obtain  by 
multiplication  of  Eq.  (51)  with  n*(k,  t-x),  after  ensemble  averaging, 

au/ax  +k2u(k,x)  +  iSk'k'<n*(k,  t-x)n(k’,t)n(k-k',t)>  =0.  (53) 

Linking  the  complex  function  U(k,x)  to  a  third-order  correlation  function,  this  is  the  first 
equation  of  an  infinite  chain  of  equations  connecting  the  order  to  the  order  N-t-1  correlation 
function.  Applying  the  operator  -a/ax  +k^  to  Eq.  (53),  and  using  the  complex  conjugate  of  Eq. 
(51)  to  define  the  action  of  this  operator  on  the  first  factor  inside  the  averaging  brackets  in  Eq. 
(53),  we  obtain 

-a2u/ax2  +  k^U  =  -ilk't<'<  [iSk"k"n*(k",t-x)n*(k-k",t-x)]n(k',t)n(k-k',t)> 

=  Zk',k"k'k''[U(k",0)U(k',0)5k.O  +U(k-,x)U(k-k',x)5k"k-  +U(k',x)U(k-k',x)5k",k-k']- 

(54) 

The  last  form  was  obtained  by  approximating  the  fourth-order  correlation  function  with  its 
expression  in  terms  of  the  second-order  correlation  applicable  for  Gaussian  processes: 

<ABCD>  =  <AB><CD>  -t-  <AC><BD>  -t-  CADXBO.  (55) 

This  approximation  was  used  by  Heisenberg  [20]  in  his  turbulence  theory,  and  was  successfully 
verified  in  its  practical  applicability  by  Uberoi  [21]  and  Batchelor  [22]. 

In  the  limit  L^oo,  setting  U(k,x)L/2x  =  u(kx),  we  write  for  k^^O  Eq.  (54)  in  the  form 

-d^u/di'^  +  k'^u  =  |[k'2u(k',x)u(k-k’,  x)  -►  k'(k-k’)u(k’,x)u(k-k',x)]dk' 


or 

-92u(k,x)/ax2  +  k4u(k,x)  =  kjk'u(k‘,x)u(k-k',  x)dk’.  (56) 

This  is  our  fundamental  dynamical  equation  of  traffic  turbulence.  All  unspecified  limits  on 
integrals  are  from  -oo  tooo,  as  we  mentioned  earlier.  From  the  symmetry  A(^,x)  =A(-^,-x) 


=A  (^,x)  we  see  that  a  physically  acceptable  solution  of  this  equation  must  satisfy  the 
conditions 


u(k,T)  =  u(-k,-x)  =  u  (-k,-x)  =  u  (k,-x). 


(57) 


A  solution  of  the  form 


u(k,x)  =  V(k)e-'mkx_  ^ith  V(k)  =  V(-k)  =  V*(k) 


(58) 


where  m  is  a  real  constant,  substituted  into  Eq.  (56),  yields 


k(m2  +  k2)V(k)  =  Jk’ V(k')  V(k-k’)dk’. 


(59 


In  the  low  wave  number  region  k«1  (i.e.,  k«vo/D)  we  neglect  the  k2  term  which  arises  from 
diffusion,  and  get 


m2kV(k)  =  jk'V(k')V(k-k')dk'. 


(60) 


A  solution  V(k)  =  QlklS'l,  with  arbitrarily  small  £>0,  yields  the  value  of  C,  independent  of  k 
only  in  the  limit  of  e-)0: 


OO  oo 

m2/C  =  Jk’|k'|e-1 11 -kVkie-IdkVk  =|k|£  jK|K|e-l|1 -kIE''' dK 

-OO  .OO 

2 

=  |k|^  J|K|fi1-K|f-''dK  -  2/e  =  I  =  const,  (6  1) 

0 

where  we  used  the  substitution  k'/'k=K.  The  divergence  at  e=0  disappears  when  we  return  to  a 
finite  value  of  the  maximal  road  length  L,  which  corresponds  to  a  minimal  k  (or  k')  value 
ko=27:/L,  and  transforms  the  integrals  back  into  sums. 

Eq.  (61)  establishes  a  proportionality  between  the  level  of  the  turbulence, described  by 
C,  and  the  magnitude  of  the  small  parameter  eas  in  earlier  turbulence  calculations  [14-17]. 
The  same  fundamental  feature  is  more  clearly  expressed,  without  the  pseudo-singularities 
present  here,  in  the  quantum  1/f  theory  [1-13],  where  e  becomes  the  infrared  exponent  aA 
known  from  quantum  electrodynamics.  Indeed,  aA  is  present  there,  just  as  we  see  it  come  in 


here,  both  as  a  factor  in  front  of  the  final  result,  and  as  a  small  defect  in  the  exponent  of  the 
frequency.  We  also  mention  that  the  apparent  absolute  determination  of  C  by  Eq.  (61)  is  just  an 
artifact  which  reminds  us  that  we  omitted  a  source  term  in  Eqs.  (47),  (49-51),  (53-54), 
(56),  which  comes  in  as  a  5  function  of  time,  h(k)5(t)  on  the  r.h.s.  of  Eq.  (53),  only  if  we 
assume  the  excitation  in  Eq.  (47)  to  depend  on  time  like  white  noise.  Once  that  source  term  is 
written  explicitly,  it  will  allow  us  to  determine  C  as  this  was  done  earlier  [14]  for  a  different 
equation,  and  it  does  not  affect  our  equations  for  0. 

We  can  rewrite  Eq.  (52)  in  the  form 

A(^,x)  ^  <n*(x,t)n(x+4,t+x)>  =  IkU(k,x)e'*<^  =  Jdk  u(k,x)e'k^ 

=  c|e'*<(^-n^'^)|k|s-''dk  =  |^-mx|-e(m2/|)Jei’^|K|e-'*dK  =  m2|^-mx|-e.  (62) 

Here  we  have  used  the  substitution  ktq-mx|=K,  and  we  have  taken  the  limit  e->0  in  the  integral 
and  in  I,  noting  that  they  exhibit  the  same  divergence,  and  therefore  can  be  simplified  in  the 
limit.  The  cancellation  of  the  divergences  in  the  expression  of  the  autocorrelation  function 
shows  that  this  important  function  is  finite  even  in  the  continuum  limit.  It  is  only  the  spectral 
distribution  of  turbulence  which  exhibits  a  singularity  at  low  wave  numbers  and  frequencies. 

Since  only  the  limit  e=0  satisfies  the  fundamental  dynamical  equation  of  traffic 
turbulence,  Eq.  (62)  means  that  the  autocorrelation  function  is  a  constant.  It  also  indicates  that 
the  constant  has  to  be  interpreted  as  the  limit  of  a  slowly  decreasing  power  law  which  depends 
only  on  |^-x|;  therefore  the  spectrum  in  wave  numbers  is  the  same  as  the  spectrum  in 
frequencies.  The  Fourier  transform  with  respect  to  x  is  the  spectral  density 

Sn'(q.w)  =J<n*(x,t)n(x+^,tt-x)>e'“^dx  =  clj(dk/ik|'' -e)e'‘^(^-n^x)+iwxdt.  (63) 

To  leave  the  moving  (primed)  frame  and  return  to  the  system  of  reference  at  rest,  we  replace  co 
by  co-kvo  in  the  last  integrand,  actually  just  by  w-k  in  our  (vo=1)  normalization: 

Sn(^.co)  =|<n*(x,t)n(x+4,t+x)>e'(‘^~*^''Hxcjt 

=  cJj(dk/k1-e)e'‘^("-'^'^)+'(w-k)xdx.  (64) 

The  spectral  density  of  concentration  fluctuations  in  a  given  point  is  obtained  by  setting  ^=0. 

The  integration  with  respect  to  x  yields  a  delta  function  27r5[a)-(m+1)k],  and  we  finally  obtain 
the  spectrum 


Sn(0,co)  =  C[(m+1)/co] 


(65) 
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where  again  the  constant  will  turn  out  to  be  finite  when  we  coarse  grain  the  integrals.  As 
mentioned  above,  the  constant  will  actually  be  determined  by  the  excitation  term  omitted  in  Eq. 
(56). 

According  to  Eq.  (2)  we  can  write  J(k,co)  =  (co/k)n(k,a))  in  the  system  at  rest. 
Therefore,  for  J  we  include  a  factor  ((o/k)^  into  the  integrand  of  Eq.  (64),  and  we  get  again  the 
spectrum 

Sj(^,w)=l<J*(x,t)J(x+^,t+T)>e'(w-kVo)T:cit 

=  clJ((d2dk/k3-E)eik(^-mt)+i(m-k)xdx: 

Sj(0,co)  =  C(m+1)2  [(m+1  )/co]'' =  Const/co. 

Due  to  the  neglect  of  the  k^  term  in  Eq.  (59),  the  universal  1/f  spectrum  will  be 
limited  towards  high  frequencies  by  axvo^/D,  and  the  1/k-spectrum  by  k<vo/D.  This 
limitation  occurs  because  only  Eq.  (60)  satisfies  the  condition  of  being  free  of  any 
characteristic  scale,  thus  exhibiting  our  sliding  scale  invariance  [23],  while  Eq.  (59)  does  not 
have  this  property.  We  have  thus  constructed  a  statistical  dynamic  theory  of  traffic  turbulence, 
proving  analytically  Musha's  earlier  result,  without  any  pretension  of  mathematical  rigor. 
Traffic  turbulence  arising  from  instabilities  of  the  laminar  traffic  flow  can  be  considered  as  a 
form  of  classical  fluid-dynamical  chaos. 

III. 4  Quantum  1/f  Noise:  Application  of  the  Sufficient 
1/f  Criterion  to  QED  Quantum  1/f  Chaos 

The  nonlinearity  causing  the  1/f  spectrum  of  turbulence  in  both  semiconductors  and 
metals  is  caused  by  the  reaction  of  the  fieid  generated  by  charged  particles  and  their  currents 
back  on  themselves.  The  same  nonlinearity  is  present  in  quantum  electrodynamics  (QED), 
where  it  causes  the  infrared  divergence,  the  infrared  radiative  corrections  for  cross  sections 
and  process  rates,  and  the  quantum  1/f  effect.  We  shall  prove  this  on  the  basis  of  the  sufficient 
criterion  for  1/f  spectral  density  in  chaotic  systems,  derived  in  the  previous  annual  report. 

Consider  a  beam  of  charged  particles  propagating  in  a  well-defined  direction  which  we 
shall  call  the  x  direction,  so  that  the  one-dimensional  Schrodinger  equation  describes  the 
longitudinal  fluctuations  in  the  concentration  of  particles.  Considering  the  non-relativistic  case 


(66) 

(67) 
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which  is  encountered  in  most  quantum  1/f  noise  applications,  we  write  in  second  quantization 
the  equation  of  motion  for  the  Heisenberg  field  operators  \\f  of  the  in  the  form 

iha\|//3t  =  (1/2m)[-ihV  -(e/c)  A]2v};,  (68) 

With  the  non-relativistic  form  J  =  -ih\}f’Vy/m  +  hermitic  conjugate,  and  with 

A(x,y,z,t)  =  (h/2cmi)J^^-|*J^“^.y'*'  ^dx’  (69) 

we  obtain 

iha^/at  =  (1/2m)[-ihV  -(eh/2c^mi)  V\|/-yy\|r 

I  X  “  X  I 

At  very  low  frequencies  or  wavenumbers  the  last  term  in  rectangular  brackets  is  dominant  on 
the  r.h.s.,  leading  to 

ih3y/3t  =  (-1/2m)  [(  e  fV2c2m)  ^dx’j^y.  (71) 

I X  -  X  I 

For  X  replaced  by  Xx,  and  x'  replaced  by  Xx‘,  we  obtain 

ih3v/3t  =  (-1/2m)[(efV2c2m)J^^----jj^'^^T"^^^  =  (72) 

This  satisfies  our  homogeneity  criterion  with  p=-2.  Our  sufficient  criterion  only  requires 
homogeneity,  with  any  value  of  the  weight  p.  for  the  existence  of  a  1/f  spectrum  in  chaos. 
Therefore,  we  expect  a  1/f  spectrum  of  quantum  current-fluctuations,  i.e.,  of  cross  sections 
and  process  rates  in  physics,  as  derived  in  detail  in  Sec.  II  above.  This  is  in  agreement  with  the 
well-known,  and  experimentally  verified,  results  of  the  Quantum  1/f  Theory. 

In  conclusion,  we  realize  that,  both  in  classical  and  quantum  mechanical  nonlinear 
systems,  the  limiting  behavior  at  low  wave  numbers  is  usually  expressed  by  homogeneous 
functional  dependences,  leading  to  fundamental  1/f  spectra  on  the  basis  of  our  criterion. 
Therefore,  we  expect  a  1/f  spectrum  of  current  fluctuations,  i.e.,  of  cross  sections  and  process 
rates  in  physics.  This  is  in  agreement  with  the  well-known,  and  experimentally  verified, 
results  of  the  conventional  Quantum  1/f  Theory. 

In  conclusion,  we  realize  that,  both  in  classical  and  quantum  mechanical  nonlinear 
systems,  the  limiting  behaviour  at  low  wave  numbers  is  usually  expressed  by  homogeneous 
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functional  dependences,  (as  shown  in  Eq.  2),  leading  to  fundamental  1/f  spectra  on  the  basis  of 
our  criterion. 


III.5  ABSENCE  OF  1/F  FLUCTUATION  SPECTRA  FOR  MODE  ENERGY  AND  PHONON 
NUMBER  IN  A  NONUNEAR  CHAOTIC  CHAIN  OF  ATOMS 

Consider  a  chain  of  atoms  in  the  x  direction,  with  a  lattice  constant  b  and  displacements 
qi  from  the  equilibrium  position.  The  equations  of  motion 

md2qn/dt2  =  A[(qn  +  i -qn)-(qn-qn-1 )]  +  B[(qn  +  i -qn)^-(qn-qn-1  )2]  + 

+  C[(qn  +  l -qn)^'(qn'qn-1  )^]  (T3) 

contain  anharmonic  terms  as  long  as  B  and  C  are  different  from  zero.  With  an=qn+i-qn  and 
Pn=qn-qn-i.  neglecting  the  second-order  term,  we  obtain 

md2qn/dt2=  [(qn  +  i -qn)-(qn-qn-1  )I[A  +  B(qn+i*qn-l)  +  C(an2+anPn+Pn^)]-  (74) 

Going  over  from  finite  differences  to  a  continuum  description,  we  obtain  a  differential  equation 
for  q(x,t) 

ma2q/at2  =  b2a2q/ax2[A  +  bBaq/ax  3Cb2(aq/ax)2].  (75) 

Performing  a  fourier  transform  with  respect  to  x, 

ma2qk/at2  =  -Ab2k2qi^  -  ib2B J(k-k')k‘2qi^-qi^.l^’dk' 

+  3Cb^|dk'jdk’k'k"(k-k'-k")2qk’qk'qk-k'-k"-  (76) 

All  integrals  are  from  minus  infinity  to  infinity.  Substituting  qk=u(k,t)exp[ikbt(A/m)''^2j^ 

ma2u/at2  2ikb(A/m)i/23u/at  =  -ib^B J(k-k-)k'2uk-Uk-k'dk' 

+  3Cb4jdk’|dk"k’k"(k-k’-k")2uk-Uk-Uk-k'-k--  (77) 

our  1/f  chaos  criterion  requires  both  nonlinearity  and  homogeneity,  as  well  as  the  presence  of 
chaos  or  of  a  quasichaotic  state.  The  nonlinearity  condition  is  satisfied  unless  B=0  and  C=0, 
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while  homogeneity  requires  the  existence  of  two  numbers  p  and  9  such  that  replacing  k  by  Xk 
everywhere  except  in  the  integration  differentials,  and  replacing  t  by  X®t  leaves  the  equation 
multiplied  by  a  general  factor  XP,  i.e.,  formally  invariant.  In  our  case  we  note  that  homogeneity 
would  be  is  satisfied  with  p=2  and  9=-1  only  if  we  neglect  both  the  second-order  and  third- 
order  terms  by  setting  B=C=0,  which  would  violate  the  nonlinearity  requirement,  in 
conclusion,  our  criterion  can  not  be  satisfied. 

Nevertheless,  in  three-dimensional  piezoelectric  crystals  quantum  1/f  fluctuations  of 
the  phonon  number  are  predicted  by  the  electrodynamic  quantum  1/f  noise  theory  as  we  show  in 
detail  in  Sec.  VI,  and  have  been  observed  experimentally  in  the  Brillouin  scattering  of  light  by 
T.  Musha.  These  quantum  1/f  fluctuations  are  predicted  by  the  sufficient  criterion  as  we  have 
seen  in  Sec.  III. 4. 

IV.  PRACTICAL  APPLICATION  TO  QUANTUM  1/f  NOISE  IN  INFRARED  DETECTORS 

Quantum  1/f  Noise  is  a  fundamental  aspect  of  quantum  mechanics,  representing 
universal  fluctuations  of  physical  process  rates  R  and  cross  sections  a  given  by  the  fractional 
(or  relative)  spectral  density  S(f)  =  2aA/fN.  Therefore  it  is  present  in  the  process  rates 
generating  the  dark  current  observed  in  junction  photodetectors,  such  as  diffusion  (scattering 
cross  sections  fluctuate)  in  diffusion-limited  junctions,  and  recombination  in  the 
recombination-limited  regime.  One  is  therefore  tempted  to  expect  similar  fluctuations  in  the 
photogeneration  of  electron-hole  pairs.  However,  as  we  show  below,  the  corresponding  quantum 
1/f  coefficient  is  zero,  precluding  the  existence  of  quantum  1/f  fluctuations  in  the  photo¬ 
generation  rate.  Here  N  is  the  number  of  carriers  used  to  define  or  measure  the  process  rate  or 
cross  section  considered. 

For  an  arbitrary  process  involving  a  total  of  n  incoming  and  outgoing  charged  particles, 
the  nonrelativistic  quantum  1/f  coefficient  is  given  {24]  by 

n 

2aA  =  (4a/37cc2)  jj]  jq  ,qj(v,-v,)2  _  (78) 

where  the  summation  runs  over  the  charges  qj  and  velocities  vj  of  all  incoming  (Tii=-1)  and 
outgoing  (Tii=1)  particles  (altogether  n  of  them)  in  the  process  whose  quantum  1/f  noise  we 
want  to  find,  and  a  is  Sommerfeld's  fine-structure  constant,  e2/iic=1/137.  In  a  photoelectric 
process  a  photon  (q=0)  is  absorbed,  and  a  pair  of  oppositely  charged  particles  is  generated 
(t|=1)  with  velocities  vi  and  V2  which  are  either  zero,  or  quickly  decay  to  zero  in  a  time 
negligible  with  respect  to  the  reciprocal  frequency  at  which  we  calculate  the  quantum  1/f  noise. 


23 


Thus  in  our  case  there  are  no  incoming  charged  particles,  and  n=0+2=2.  The  aA  coefficient  of  a 
photogeneration  process  is  therefore  zero, 

aAph  =  (1,1)+(2,2)+(1,2)+(2,1)  =  0+0+  (4a/37rc2)(vi -V2)2  =  0.  (79) 

All  photogenerated  carriers  of  the  right  sign  are  collected  in  the  well  of  the  charge-coupled 
device,  although  they  may  generate  quantum  1/f  voltage  fluctuations  on  their  way.  Since 
usually  only  the  number  of  carriers  collected  at  read-out  matters,  no  quantum  1/f  noise  will  be 
observed  in  a  photoelectric  CCD  as  long  as  the  dark  current  is  negligible  with  respect  to  the 
photocurrent.  This  is  in  agreement  with  the  experiments  performed  by  Mooney  [25]  of  RADC- 
Hanscom  AFB.  The  same  considerations  apply  to  Metal-Insulator-Semiconductor  (MIS) 
photodetectors  [26]. 

V.  QUANTUM  1/F  NOISE  IS  QUANTUM  CHAOS:  FRACTAL  DIMENSION 
OF  QU/\NTUM  1/F  NOISE 

Quantum  1/f  noise  in  a  physical  quantity  such  as  a  current  j,  a  cross  section,  a  process 
rate,  or  a  kinetic  coefficient,  such  as  the  mobility,  is  represented  by  an  expression  of  the  form 

5]/<j>  =  a  Xib(k,X)|2cos(ckt+YKx)  (80) 

k;X 

where  the  sum  is  over  all  electromagnetic  modes  labeled  by  their  wave  vector  k  and  polarization 
k,  with  k=|k|.  Eq.  (80)  is  deterministic  and  the  particle-specific  random  initial  phases  are 
present  in  each  term  as  initial  conditions,  like  the  initial  phases  describing  a  turbulent  fluid 
are  present  in  each  Fourier  component  of  the  perfectly  deterministic  velocity  field.  The 
difference  between  the  mentioned  classical  chaos  (turbulence  )  case  and  the  quantum  1/f  chaos 
we  are  introducing  here  becomes  evident  when  we  recall  that  J  is  in  fact  a  probability  current 
density!  Can  we  calculate  a  spectral  density  of  the  fluctuations  5]  of  this  probability  current 
and  claim  that  it  represents  the  expectation  value  of  the  spectral  density  of  the  observed 
quantum  1/f  fluctuations?  Our  rigorous  derivation  in  second  quantization  [23,  28-29]  tells  us 
we  can,  provided  we  divide  the  single-particle  result  by  N  in  the  case  of  bosons  and  N-1  for 
fermions.  Here  N  is  the  number  of  scattered  particles  used  to  define  what  we  call  scattered 
current  j  or  scattering  cross  section  o.  Encouraged  by  this  result,  we  have  undertaken  an  effort 
to  determine  the  fractal  dimension  of  quantum  1/f  noise  given  by  Eq.  (80)  and  restricted  to  an 
observed  frequency  interval  from  cko=  eo=2;tfo  to  a=2jiF,  where  eo=27:/T  is  determined  by  the 

duration  T  of  the  1/f  nose  measurement.  Our  objective  is  to  determine  the  fractal  dimension  of 


quantum  1/f  noise  and  to  compare  It  with  measurements  of  the  fractal  dimension  of  physical  1/f 
noise  in  HgCdTe  MWIR  photodiodes  [30],  Calculating  the  spectral  density,  replacing  the 
summation  In  Eq.  (80)  by  an  Integral,  and  going  back  to  the  original  fluctuations,  we  obtain  an 
equivalent  representation  of  the  fluctuations  (23,28-29) 

A 

5J/<J>=i  2(aA)’'2  jQos(et+Ye)deNe  (81) 

Co 


where  a=1/137  Is  Sommerfeld's  fine  structure  constant  and  aA  ■  (2a/3n)(Av/c)2  is  the 
bremsstrahlung  coefficient,  or  Infrared  exponent,  of  the  process  which  generates  the  current  J, 
Av  being  the  velocity  change  of  the  charged  particles  In  the  process.  Repeating  the  fluctuations 
periodically  outside  the  observational  Interval  T,  we  can  represent  Eq,  (81)  through  an 
equivalent  Fourier  series 


M 

5j/<j>«  2aA  ]^cos(neot+Yn)/(neo)’'^  ■  x(t)  (8  2) 

0=1 

which  Is  similar  to  Eq.  (80),  with  the  terms  of  the  same  frequency  grouped  together,  and  the 
terms  lumped  together  In  harmonics  with  random  phases  Vn.  with  M<=A/eo.  To  this  expression 

we  have  applied  the  Takens-Grassberger-Procaccia  analysis  [31],  by  creating  a  time  series 
xi=x(ti)  with  l<kN.  We  consider  this  series  as  a  one-dimensional  sample  of  the  quantum  1/f 
process.  We  calculate  the  correlation  function  Ci(*^)(r)  which  Is  N‘2  times  the  number  of  data 
pairs  (XI, xj)  separated  by  a  distance  less  than  r.  Next  we  form  groups  of  d  consecutive  data 
(xi,xj+i,.,.xi+d.i).  Considering  them  as  vectors  In  a  d-dimensional  Euclidean  space,  we  again 
calculate 

N 

CdC^^fr)  -  N-2  5^0[r-llxi-xj||].  (83) 

1  ;j*  1 

where  0(x)  is  the  step  function  (zero  for  negative,  1  for  positive,  and  1/2  for  null  arguments) 
Finally  we  plot  the  curves  logioCd^^^r)  as  a  function  of  r,  with  d=0, 1,  2,.„  as  a  parameter. 
We  notice  (Fig.  1  for  M=*1  and  Fig.  2  for  M=2  and  3  where  for  M=:2  the  numbers  d=1 ,2.3,4, 
have  been  replaced  by  a,b,c,d.)  that  the  slope  of  the  curves  Is  Just  d,  and  Increases  as  d  Is 
increased  from  curve  to  curve  up  to  M,  the  number  of  terms  In  Eq.  (82).  We  graphed  the  case 
M=1  (a  single  term  in  Eq.  (82)  In  Fig.  1  and  M=2  and  3  In 
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Fig.  2.  After  d=M,  the  slope  does  no  longer  increase.  This  maximal  slope  gives  the  fractal 
dimension  of  the  quantum  1/f  process  \with  M  terms.  Our  conclusion  at  this  point  is  that 
quantum  1/f  noise  in  a  frequency  interval  corresponding  to  M  terms  in  Eq.  (83)  is  chaotic  with 
a  fractal  dimension  M.  Had  it  been  stochastic  rather  than  chaotic,  the  slope  of  the  curves  would 
have  kept  increasing  indefinitely. 

Experimentally,  Fote,  Kohn,  Fletcher  and  McDonough  found  a  fractal  dimension  of  dM=10 
for  1/f  noise  measured  in  HgCdTe  MW  infrared  detectors  in  the  interval  from  fo=10"2Hz  to 
F=10  Hz.  Noticing  that  F/fo=1 0^=210,  we  expect  M=10,  which  yields  also  dM=10.  This  nice 
agreement  between  theory  and  experiment  indicates  that  the  measured  1/f  noise  is  a  form  of 
(deterministic)  chaos,  whereas  the  other  models  which  competed  with  the  quantum  1/f  theory 
incorrectly  described  it  as  a  stochastic  phenomenon.  Furthermore,  this  agreement  represents 
an  independent  verification  of  the  quantum  1/f  theory,  independent  from  tests  based  on  the 
predicted  magnitude  and  spectral  dependence  of  1/f  noise. 


Vi.  QU/WTUM  1/f  FUCKER  OF  FREQUENCY  IN  QUARTZ  RESONATORSTHEORY  AND  EXPERIMENT 


VI. 1.  Introduction 


Flicker  of  frequency  noise  is  an  important  characteristic  of  quartz  resonators  which 
limits  their  stability  and  determines  their  utility  for  most  applications.  The  1/f  contribution 
to  frequency  stability  is  best  obtained  by  observing  the  stability  over  a  range  of  measurement 
times  in  the  time  domain,  or  over  a  range  of  frequencies  in  the  frequency  domain.  From  the 
extended  data  one  can  fit  a  flicker  of  frequency  model  to  the  data  that  excludes  the  contributions 
from  random  walk  frequency  modulation  and  from  the  drift  present  in  many  resonators  and 
oscillators  (Fig.  3). 

Qne  of  the  first  who  systematically  studied  the  1/f  noise  as  a  function  of  geometry, 
temperature  and  Q-factors  was  Gagnepain  [32].  He  noticed  empirically  that  the  1/f  part  of  the 
spectral  density  of  fractional  frequency  fluctuations,  Sy(f),  varied  as  Q''^  for  resonators 
between  1  and  25  MHz.  As  the  temperature  changes,  the  Q-factor  of  a  resonator  changes.  This 
allows  us  to  exclude  the  effect  of  many  other  factors.  Additional  work  by  Parker,  however, 
showed  that  the  data  from  both  bulk  acoustic  wave  (BAW)  and  surface  acoustic  wave  (SAW) 
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devices  could  be  roughly  fit  to  the  the  same  model  [33],  if  one  assumes  a  Q'^  dependence  for  the 
phase  noise  S,),  insted  of  Sy.  The  fit  is  shown  on  Fig.  4. 

From  a  theoretical  point  of  view,  fundamental  work  by  this  author  [34],  inspired  by  the 
quantum  1/f  theory  [2],  has  derived  the  depenence  of  Sy  for  all  resonators  from  first 
principles  already  1979.  Work  on  many  systems  other  than  quartz  has  yielded  very  good 
quantitative  agreement  between  theory  and  experimental  data  for  quantum  1/f  noise  [35].  This 
paper  refines  the  previous  theoretical  work  on  1/f  frequency  noise  in  quartz  to  suggest  a  better 
framework  for  predicting  the  level  of  1/f  frequency  noise  in  quartz  resonators  over  a  wide 
range  of  frequencies  and  Q-factors. 


Vi.  2.  Theory  of  1/f  Frequency  Noise  in  Quartz  Resonators 

According  to  the  general  quantum  1/f  formula  [2],  r“2Sp(f)=2a/Vf  with  a=e2/hc 
=1/137  and  A=2(AJ/ec)2/37i  is  the  quantum  1/f  effect  in  any  physical  process  rate  r.  Setting 
J=dP/dt=P  where  P  is  the  vector  of  the  dipole  moment  of  the  quartz  crystal,  we  obtain  for  the 
fluctuations  in  the  rate  r  of  phonon  removal  from  the  main  resonator  oscillation  mode  (by 
scattering  on  a  phonon  from  any  other  mode  of  average  frequency  <©>)  of  the  crystal,  (or  via  a 
two-phonon-process  at  a  crystall  defect  or  impurity,  involving  a  phonon  of  average  frequency 
<co'>)  the  spectral  density 

Sr(f)  =  r24a(AP)2/37te2c2,  (84) 

where  (aP)^  is  the  square  of  the  dipole  moment  rate  change  associated  with  the  process  causing 
the  removal  of  a  phonon  from  the  main  oscillator  mode.  To  calculate  it,  we  write  the  energy  W 
of  the  interacting  resonator  mode  <ci)>  in  the  form 

W  =  nh<a)>  =  2(Nm/2)(dx/dt)2  =  (Nm/e2)(e  dx/dt)2  =  (m/Ne2)E2(P)2  ;  (85) 

The  factor  two  includes  the  potential  energy  contribution.  Here  m  Is  the  reduced  mass  of 
the  elementary  oscillating  dipoles,  e  their  charge,  e  a  polarization  costant,  and  N  their  number 
in  the  quartz  crystal.  Applying  a  variation  An=1  we  get 
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An/n  =  2|aP|/|P|,  or  AP=P/2n. 

Solving  Eq.  (85)  for  P  and  substituting,  we  obtain 


(86) 


|AP|  =  (NfKco>/n)l/2(e/2e)  (87) 

Substituting  AP  into  Eq.  (84),  we  get 

r'2Sr(f)  =  Nati<co>/3n7cmc2f£2  =  A/f  .  (88) 

This  result  is  applicable  to  the  fluctuations  in  the  loss  rate  r  of  the  quartz. 

The  corresponding  resonance  frequency  fluctuations  of  the  quartz  resonator  are  given  by2 

co-2So)(f)  =  (1/4Q)(A/f)  =  Nah<u>/12njtmc2f£2Q4^  (89) 

where  Q  is  the  quality  factor  of  the  single-mode  quartz  resonator  considered,  and  <co>  is  not  the 
circular  frequency  of  the  main  resonator  mode,  coq,  but  rather  the  practically  constant 
frequency  of  the  average  interacting  phonon,  considering  both  three-phonon  and  two-phonon 
processes.  The  corresponding  aP  in  the  main  resonator  mode  has  to  be  also  included  in 
principle,  but  is  negligible  because  of  the  very  large  number  of  phonons  present  in  the  main 
resonator  mode. 

Eq.  (89)  can  be  written  in  the  form 

S(f)  =  pV/fQ4,  (90) 

where,  with  an  intermediary  value  <(d>=10®/s,  with  n=kT/fi<a)>,  T=300K  and  kT=4 

p  =(N/V)ah<Q)>/12n7:e2mc2  =1022(1/137)(1  o8)2/i  2kT7i1  0-27  g  1020  =-|  (91) 

The  form  of  Eq.  (30)  shows  that  the  level  of  1/f  frequency  noise  depends  not  only  as  Q"'^ 
as  previously  proposed,  but  also  on  the  oscillation  frequency  or  the  volume  of  the  active  region. 
This  model  qualitatively  fits  the  data  of  Gagnepain  [32,42  ]  where  he  varied  the  Q-factor  with 
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temperature  in  the  same  resonator  (but  not  frequency  or  volume). 

The  model  also  provides  the  basis  for  predicting  how  to  improve  the  1/f  level  of 
resonators,  beyond  just  improving  the  Q-factor,  which  has  been  known  for  many  years.  Since 
the  level  depends  on  active  volume,  one  should  use  the  lowest  overtone  and  smallest  diameter 
consistent  with  other  circuit  parameters.  The  following  experimental  discussion  has  been 
generated  by  F.  Walls  (NIST  Boulder). 


VI.  3.  Experimental  Measurements  and  Analysis  of 
1/f  Noise  in  Quartz  Resonators 

The  level  of  1/f  frequency  noise  in  quartz  resonators  has  been  measured  using  phase 
bridges  and  complete  oscillators  [32,33,36-41].  Unfortunately  much  of  the  data  in  the 
literature  is  unusable  for  modeling  because  the  unloaded  Q-factor  is  unknown.  (Our  case  is 
even  more  restrictive  because  we  also  need  to  know  the  electrode  size).  The  phase  bridge 
approach  has  the  advantage  that  the  unloaded  Q-factor  can  be  easily  measured  and  the  noise  in 
the  measurement  electronics  can  be  evaluated  independent  of  the  resonator.  If  resonator  pairs 
are  used,  the  noise  of  the  driving  source  can  generally  be  neglected  and  the  pair  can  operate  at 
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Function  of  Loaded  Q. 


virtually  any  frequency  [39],  The  oscillator  approach  makes  it  possible  to  compare  many 
different  resonators  one  at  a  time.  The  noise  of  individual  oscillators  can  be  derived  by 
measuring  the  phase  noise  between  3  oscillators  [42]. 

Figure  5  taken  from  [33]  is  one  of  several  studies  showing  that  the  1/f  level  is 
virtually  independent  of  the  loaded  Q-factor.  This  is  in  complete  agreement  with  the  theoretical 
model.  In  practical  oscillators  there  is  a  dependence  on  loaded  Q-factor  only  when  the  phase 
noise  of  the  sustaining  electronics  contributes  to  the  overall  noise  level. 

We  have  analyzed  1/f  frequency  noise  as  a  function  of  unloaded  Q,  volume  under  the 
electrodes,  and  frequency.  For  a  given  resonator  geometry  and  manufacturer,  we  have  taken  the 
best  values  of  Sy(f)  reported  in  an  attempt  to  remove  the  effects  of  poor  crystals  or 
electronics.  In  Fig.  6  we  have  taken  all  of  the  precision  data  available  with  unloaded  Q-factor, 
electrode  volume,  and  frequency  stability  and  plotted  them  according  to  the  three  models.  Except 
for  the  2.5  MHz  resonator  where  Qvq  =  0.95x1  the  Qvq  product  for  all  resonators  is  near 

1.2x10''2  (this  is  close  to  the  material  limit  for  AT  and  SC  cut  resonators).  The  curve  labeled 
Ky  shows  the  fit  of  the  data  to  the  model  [32] 

Sy  (f)  =Ky  /f(3  X  10'^  /Q^  ).  Ky  vories  about  a  factor  of  500  for  Q-factors  between  10^  and 

3.8x10®  (resonator  frequencies  between  2.5  and  100  MHz). 

The  curve  labeled  K(j  shows  the  fit  of  the  same  resonator  data  to  the  model  [33] 

S  (f)  =  Ko/f®  (3  X  10’°  /Q'‘  ).  Kq  varies  about  a  factor  of  10  for  the  same  range  in  Q-factor. 
Curves  Pe  s^d  pb  show  the  fit  of  the  same  resonator  data  to  the  model  Sy(f)  =  p/f  (Vol/Q'*  ), 
where  Pe  is  for  SC  and  AT  resonators  with  electrodes  plated  on  the  resonator  and  Pb  is  for  BVA- 
style  AT  and  SC  resonators  [37].  The  volume  between  the  electrodes  (in  cm^)  is  used  to 
approximate  the  volume  of  quartz  contributing  to  the  output  power.  The  B  factors  are 
remarkably  constant  for  Q-factors  from  1Q®  J2.MJLIO®. 

Figure  7  shows  the  dependence  of  Pe  on  Q-factor  for  a  number  of  electroded  resonators  of 
the  same  type  from  a  single  manufacturer  for  3  resonator  types  as  measured  by  Norton 
[38,45].  The  wide  variation  in  pg  for  the  same  style  resonator  and  Q-factor  indicates  that 
acoustic  loss  is  not  the  only  mechanism  contributing  to  the  noise  level.  The  data  for  this  graph 
was  taken  from  measurements  of  (100  s)  and  may  have  been  biased  high  by  random  walk  FM 
noise  in  some  resonators. 

VI.  4.  Discussion 

The  analysis  of  the  most  stable  quartz  resonators  indicates  that  the  1/f  frequency  noise 
level  depends  on  volume  between  the  electrodes  and  unloaded  Q-factor  in  relatively  good 
agreement  with  eq.  (30)  considering  the  fact  that  the  estimation  of  <co>  and  N  is  only 
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approximative.  The  nonadjustable  parameters  Pe  and  Pb  are  virtually  constant  versus  unloaded 
Q-factor,  which  is  in  stark  contrast  to  fitting  parameters  and  Ky.  It  is  not  surprising  to  us 
that  Pe  and  Pb  are  different  for  the  two  types  of  resonators,  since  energy  trapping  and  electrode 
stress  are  considerably  different.  Fig.  7  shows  that  there  are  other  noise  processes  besides 

acoustic  losses  that  affect  the  1/f  noise  level  in  some  resonators. 

Although  we  have  analyzed  only  the  data  for  a  few  resonators,  the  consistency  of  Pe  and  Pb 

over  a  factor  of  40  in  Q-factor  and  resonator  frequency  and  the  general  agreement  for  the 
magnitude  of  p  between  theory  and  experiment  give  us  some  confidence  that  this  new  model  can 
be  used  to  predict  the  best  performance  of  different  resonator  geometries. 

This  new  volume  model  predicts  that  a  resonator  having  smaller  electrodes  would  have  a 
lower  level  of  1/f  frequency  noise  than  another  one  with  the  same  frequency  and  Q-factor  but 
larger  diameter  electrodes.  The  decrease  in  electrode  area  would  increase  the  impedance  levels 
and  degrade  the  wide-band  noise  somewhat.  For  most  resonators  the  wideband  noise  is  dominated 
by  the  electronics  and  not  the  resonator.  The  increase  in  series  resistance,  obtained  by 
decreasing  the  electrode  area  by  a  factor  of  4,  would  probably  be  tolerable  from  the  standpoint 
of  wideband  noise  but  might  require  a  change  in  loop  gain. 

BT  resonators  are  potentially  interesting  in  that  they  offer  a  Qvq  product  approximately 

three  times  higher  than  that  of  AT  and  SC  resonators.  BT  cuts  are  roughly  as  sensitive  to 
temperature  transients  as  AT  cuts.  Therefore  to  achieve  parts  in  10'’“*  frequency  stability  with 
BT  cuts  would  require  temperature  stabilities  of  order  10'®  K/s  or  100  times  better  than  is 
required  for  SC  cut  resonators  [45]. 

Based  on  these  early  observations  it  appears  that  the  level  of  1/f  frequency  noise  in 
quartz  may  yet  be  improved  to  the  low  10  ’^  level  by  applying  one  or  more  of  the  following 
techniques;  reducing  the  electrode  area,  using  BVA  type  resonators,  going  to  lower  frequencies, 
using  BT  cut  resonators.  It  must  be  remembered  that  acceleration  induced  effects  become  more 
dominant  as  the  stability  improves. 


VII.  EXPERIMENTAL  CHECKS  ON  COLLECTOR  QUANTUM  1/f  NOISE  IN  BJT'S 

VII. 1  Introduction 

1/f  noise  in  bipolar  junction  transistors  (BJTs)  was  treated  by  van  der  Ziel  [46]-[48] 
who  applied  a  Hooge-type  approach  similar  to  Kleinpenning's  treatment  [49]  of  pn  junctions, 
and  used  experimental  data  to  determine  the  Hooge  constant  which  was  in  turn  compared  with  the 
quantum  1/f  theory.  However,  since  the  BJT  is  a  minority  carrier  device,  it  requires  the 
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application  of  the  quantum  1/f  equation  [1,2]  from  the  beginning,  for  the  correct 
interpretation  of  the  number  of  carriers  in  the  denominator  of  the  Langevin  noise  source. 

In  the  most  elementary  model  [50]  of  a  BJT,  the  collector  current  Ic  arises  from 
minority  carriers  injected  from  the  emitter  into  the  base,  which  diffuse  across  the  width  Xb  of 
the  base  and  are  then  all  swept  across  the  reverse-biased  collector  junction  by  the  built-in 
field  of  the  junction.  If  we  neglect  the  usually  small  leakage  current  of  the  collector  junction 
and  the  small  fraction  of  the  carriers  recombining  in  the  base,  we  get  for  a  n+pn  BJT: 

lc=  AqDn[noB(exp(qVBE/l<T))/XB],  (92) 

where  A  is  the  cross  sectional  area  of  the  base,  q=-e  is  the  charge  of  the  minority  carriers  in 
the  base.  Dp  their  diffusion  coefficient  in  the  base,  nB(0)=noBexp(qVBE''l<T)  is  the  electron 
concentration  at  the  limit  of  the  emitter  space  charge  region,  Vbe  is  the  applied  base  -  emitter 
voltage,  and  Xb  is  the  width  of  the  base.  The  expression  in  rectangular  brackets  is  the  electron 
concentration  gradient  calculated  with  the  boundary  condition  of  a  vanishing  electron 
concentration  at  the  limit  of  the  collector  space  charge  region.  We  assume  the  base  to  be  much 
narrower  than  the  electron 

diffusion  length  Ln=(Dn-c),  XB<<Ln.  but  sufficiently  wide  to  avoid  ballistic  electron  transport 
across  the  base.  Usually  Xb  is  a  fraction  of  a  micron. 

Quantum  1/f  fluctuations  of  the  collisional  cross  sections  of  the  electrons  in  the  base 
will  yield  fluctuations  of  the  diffusion  constant,  and  of  the  mobility  (5Dn/Dn=5u/|i): 

51c  =  Aq(5Dn)[noBexp(qVBE/kT)/XBl.  (93) 

The  corresponding  spectral  density  of  fractional  fluctuations  |■2S|c  is 

lc'2<(5lc)^>f  =  Dn'2<(5Dn)2>  =  u'2<(54)2>  =  ap/fN.  (94) 

In  the  last  step  our  quantum  1/f  equation  [1,2]  was  used,  where  N  is  the  number  of  carriers 
which  define  the  scattered  (or  diffused)  current,  leaving  the  base  and  emerging  in  the  collector, 
while  Op  =  aAp  is  the  effective  quantum  1/f  noise  coefficient,  or  Hooge  constant.  The  number  of 
electrons  N  is  thus  determined  by  the  effective  lifetime  xc  of  the  electrons,  which  will  be 
slightly  lower  than  the  lifetime  in  the  unbounded  collector  material,  due  to  the  collector  lead 
contact  processes,  and  due  to  lateral  surface  recombination.  Indeed,  we  can  write  N  =  xcic/q- 
Thus  we  finally  obtain  the  spectral  density  of  the  collector  current  fluctuations: 
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Sic  =  ctnlcC|/(f'^c). 


(95) 


in  which  xc  is  the  effective  lifetime  of  the  majority  carriers  in  the  collector.  This  expression 
is  simpler,  but  similar  to  the  expression  derived  earlier,  with  the  important  difference  that 
now  we  have  a  lifetime  of  the  carriers  in  the  denominator,  while  before  it  was  the  usually  much 
smaller  diffusion  time  Td=XB^/Dn  of  the  electrons  in  the  base.  Eq.  (95)  also  implies  that  in 
narrow-base  BJTs  of  various  base-widths  an  will  be  constant,  as  in  other  devices,  rather  than 
an/td-  In  the  following  section  we  show  that  this  expression  is  in  good  agreement  with  the 
experimental  data  in  BJTs  with  1/f  collector  noise  spectra. 

VII. 2  Comparison  of  the  Calculated  Results  with  Experimental  Data 

1.  The  effective  Quantum  1/f  noise  coefficient, 
or  Hooae  constant,  an 

We  consider  the  following  scattering  processes  [51]  in  the  calculation  of  the  quantum 
1/f  noise  coefficient  an: 

a)  Normal  collision  processes  (Impurity  scattering,  Optical  scattering  and  Acoustical 
phonon  scattering)  [10]. 

b)  Intervalley  scattering:  there  are  two  types  intervalley  processes,  i.e.,  g-processes 
which  include  Umkiapp,  and  f-processes  [10']. 

From  these  points  of  view,  we  obtain  the  current  spectral  density  in  the  form  of  Eq. 

(95); 


S|c(f)  =  an[lcC|/(f'tc)l  =  aAn[lcQ^(ftc)] 

For  the  case  a),  the  normal  collision  processes: 

an  =aAn  =a4Av2/(3”C^)  =4a(kT/mo:tC^)  =4.69x10'^°,  (96) 

where  we  used  a  =  \ioCe^/2h  «  1/137  (the  fine  structure  constant),  c  =  3x10®m/sec,  k  = 
1.38x1  0-23 J/K,  J  ^  300K,  and  mo  =  g.lxIO-^^kg. 

For  the  case  b),  the  intervalley  scattering  +  umkiapp  scattering  1/f  noise  (g- 
processes): 

an  =  uAn  =  a4Av2/(37tc2)  =  a4(fiAk/m)2/3TCc2  =5.86x10"^,  (97) 
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where  av  =  Ap/m  =  Mk/m,  |Ak|  =  Q.Q{2nla),  and  a  =  5.4A  for  Silicon  [51].  The  g-proceses 
include  Umklapp  back  to  the  original  Brillouin  zone,  and  the  conduction  effective  mass  m  = 
0.26mo  must  be  used  [lO’j. 

Comparing  the  an  =  5.86  x  IQ-"^  with  reference  [10]  Fig.  8  where  the  result  was  from 
the  exact  calculation,  we  find  out  that  they  are  pretty  close.  We  would  like  to  point  out  that  an 
=  aintervalley  which  is  only  an  approximation  due  to  the  high  aintervalley  comparing  with 
^impurity  3nd  aacoustic. 

2.  The  experimental  data  compared  to  the  Hooae  parameters  ann 
i.  In  n+-p-n  bipolar  transistors: 

S^(f)/|i2  =  an/fN.  (98) 

This  can  be  written  in  terms  of  the  diffusion  constant  ,  D  =  kT)x/q 

SD(f)/D2  =  ttH/fN  (99) 


which  yields  [52] 


S|c(f) 


gHn  qlcDn  I  ,N(0) 
f  wb^  ^N(wb)'’ 


(100) 


where  aH  is  the  Hooge  parameter  for  electrons,  f  is  the  frequency,  Dn  is  the  diffusion  constant 
for  electrons,  wb  is  the  width  of  the  transistor  base  region,  N(0)  is  the  electron  concentration 
for  unit  length  at  the  emitter  side  of  the  base,  and  N(wb)  is  the  electron  concentration  for  unit 
length  at  the  collector  side  of  the  base. 

If  we  introduce  the  ratio  [53] 


N(0)  ^Vn+Pn/WB 
N(wb)  ■  D  n/WB 


(101) 


where  Vn  is  the  saturation  velocity  of  the  electrons  in  the  base  region,  and  the  diffusion  time 

Xdn  =  WB2/2Dn  =  1/27ifT  (102) 

where  fj  is  the  upper  cut  off  frequency  of  the  BJT,  then 


(103) 


S|c(f)  = 


ttHn  QIc 

f  2Xtjn 


Vn  +  (Dn/WB) 

D„/w,  1  ■ 


aHn 

f 


qlcicfTln[ 


Vn+Dn/w 

Dh/Wb 


B 


•] 


ii.  In  p+-n-p  bipolar  transistors: 


S|c(f)  = 


*^*cDpi  ,  P  ( Q.).. 

f  Wb^  P(wb) 


where 


P(0I_1 
P(Wb)  ■ 


+  and  Tdp  =  Wb2/2Dp, 


(104) 

(105) 


and  where  vcp  (=10^  cm/s  in  Si)  is  the  saturation  velocity  of  holes  in  the  base.  For  details 
see[53],  and  the  following  comparisons  in  table  I,  where  a  refers  to  the  normal  scattering 
quantum  1/f  calculation  and  b  to  the  more  likely  case  of  g-type  intervalley-umklapp  scattering. 
It  is  more  likely  to  have  this  case  since  it  has  a  much  larger  quantum  1/f  effect,  and  will  mask 
the  smaller  contribution  calculated  for  case  a.  Also,  the  experimentally  noticed  strong  increase 
in  1/f  noise  if  the  transistors  are  cut  from  single-crystals  so  that  the  current  flows  along  an 
[100]  -  like  direction  confirms  this  and  can  not  be  explained  without  the  quantum  1/f  theory. 
The  experimental  values  in  Table  I  are  much  closer  to  the  values  calculated  for  the  case  b 
replacing  ipd  with  the  lifetime  in  the  collector,  tc.  This  corresponds  to  the  inclusion  of  a 
corrective  factor  xnd/tc  in  the  derivation  of  the  collector  quantum  1/f  noise  spectral  density. 
However,  a  proportionality  of  Sic  to  fj  is  noticed  in  many  transistors,  contrary  to  our 
suggestion  of  replacing  tnd  with  xc-  The  present  formula  containing  xc  is  applicable  to 
ultrasmall  BJTs  with  a  very  narrow  base  region  for  which  most  of  the  life  time  of  carriers 
diffusing  from  the  emitter  is  spent  in  the  collector.  This  subject  requires  further  theoretical 
and  experimental  study  to  determine  the  exact  limits  of  applicability  of  the  new  formula. 
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3.76E-25 

4.69E-22 


Vlll.  QUANTUM  1/f  MOBILITY  FLUCTUATIONS  IN  SEMICONDUCTORS, 

CALCULATED  WITH  THE  QUANTUM  1/f  CROSS-CORRELATION  FORMULA 

Vlll.  1.  introduction 

Together  with  the  graduate  student  Thomas  H.  Chung,  the  author  has  performed  an 
analytical  calculation  of  mobility  fluctuations  in  silicon  and  gallium  arsenide,  using  the  new 
quantum  1/f  cross-correlations  formula  derived  by  the  author  in  the  previous  AFOSR  Grant 
period,  (1987)  and  included  in  the  July  1989  Final  Technical  Report.  This  same  new  form  of 
the  quantum  1/f  cross  correlations  was  rederived  with  a  different  method  by  Van  Vliet  in 
1989.  It  differs  from  the  old  form  used  in  the  1985  calculation  of  Kousik  et.  al.  by  a  correction 
which  is  zero  when  the  momentum  changes  of  the  two  current  carriers  involved  in  the  cross 
correlation  are  identical,  but  increases  to  finite  values  when  the  momentum  differences  caused 
by  the  scattering  process  are  different.  The  correction  is  proportional  to  the  squared  difference 
of  the  two  momentum  changes.  We  have  repeated  all  calculations  in  the  original  paper  by  Kousik 
et.al.  [10,  10'],  obtaining  both  for  impurity  scattering  and  for  the  various  types  of  phonon 
scattering  new  analytical  expressions  which  show  a  considerable  increase  of  the  final  quantum 
1/f  noise.  The  results  obtained  are  applicable  both  to  direct  and  indirect  bandgap 
semiconductors. This  calculation  is  of  major  importance  for  the  1/f  noise-related  optimization 
both  of  the  two  types  of  materials,  and  of  the  many  devices  constructed  with  them  for  military 
and  civilian  applications  in  the  electronic  and  opto-electronic  industry. 

The  new  cross-correlation  formula  gives  the  cross-spectral  density  which  describes  the 
way  in  which  symultaneous  quantum  1/f  scattering  rate  fluctuations  AW  obseved  in  the 
direction  of  the  outgoing  scattered  wave-vector  K‘  are  correlated  with  those  in  the  K"  direction, 
when  the  two  corresponding  incoming  current  carriers  have  the  wave  vectors  Ki  and  Kz: 

S^w(Ki,K';K2.K";f) 

=  (2a/3;if)(h/m'c)2WKi,K'WK2.K-[(K'-Ki)2  +  (K"-K2)2]5Ki.K2-  (106) 

The  form  conjectured  by  us  earlier  had  2(K’-Ki)(K"-K2)  in  place  of  the  rectangular  bracket. 
The  difference  between  the  rectangular  bracket  and  2(K'-Ki )(K"-K2)  is  the  perfect  square 
[(K'-Ki )-(K"-K2)]^.  Therefore  we  expect  the  new  results  to  be  always  larger  than  the 
results  obtained  on  the  basis  of  the  previously  conjectured  form. 

Vlll.  2.  Impurity  Scattering 

For  impurity  scattering  of  electrons  in  solids,  fluctuations  Ax  of  the  collision  times  x 
will  cause  mobility  fluctuations 
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A^band(t)  =  [e/m*«v2»]lKVK^AT:(t)nK. 


(107) 


where  «v2»  is  both  the  average  over  all  states  of  wave-vectors  K,  with  occupation  numbers 
nK,  in  the  conduction  band,  and  the  thermal  equilibrium  average  of  the  quadratic  carrier 
velocities.  With  the  help  of  the  relation 

1/t(K)  =  (V/8ti:0)J(1  -  cos07cos0} WK.K'd^K',  (108) 

the  mobility  fluctuations  are  reduced  to  fluctuations  of  the  elementary  scattering  rates  Wk,K'. 
governed  by  Eq.  (106).  Here  V  is  the  volume  of  the  normalization  box  which  disappears  in  the 
final  result,  0  and  0’  respectively  the  angles  K  and  K'  form  with  the  direction  of  the  applied 
field.  One  finally  obtains  after  tedious  multiple  integrations 

ti'^SA4(f)  =  [256::aK-e-*hi2/3m'8z4e8Ni2](l/f)iKKl  °[ln(1-ha2)-a2/(1  +a2)]-3 

[(2a2  +  a"^)/(1 +a2)-2ln(1 +a2)]F(EK)[2KVK^T:(K)F(EK)]'^.  (1  09) 

where  a=2K/K,  K2=e2n(T)/£kBT.  n(T)  is  the  electron  concentration,  F(EK)=exp(EF-EK)  for 
non-degenerate  semiconductors,  Nj  the  concentration  of  impurities  of  charge  Ze  and  e  the 
dielectric  constant.  The  corresponding  partial  Hooge  parameter  for  impurity  scattering  is  thus 

oo 

tti  =  [4^7:aKh5Nc/3m'7/2(kgT)3'2c2]  jdxxl  1/2e-!< 

0 

[ln(bx-h1)-bx/(bx  +  1)]-0((2bx-^b2x2)/(bx-^1)-2ln(bx-Kl)] 

OO 

{jdxx2e'’‘[ln{bx-t-l)-bx/(bx-t-l)]'’}’2.  (110) 

0 

This  result  is  graphed  below  for  three  different  values  of  the  donor  concentration  Nu  (=1021, 
1023_  and  1024)  3njj  j5  compared  with  old  results  (marked  with  an  initial  letter  o)obtained  by 
simply  recalculating  the  old  analytical  expression  (10). 


45 


Alpha  Impurity 


Temp.  (K) 


As  expected,  the  new  cross  correlation  formula  leads  to  higher  a,  values  than  the  previously 
conjectured  expression.  This  was  mentioned  in  connection  with  Eq.  (106)  above. 

Viil.  3.  Acoustic  Electron-Phonon  Scattering 

In  this  case  the  calculation  is  similar,  and  leads  to  the  result 

oo 

oac  =  [32;iaNcm*C7h3/3c2kBT)4]{(1/R2)  |dxx-'» 

1 

[(x-1)7/7  +  {R+1){x-1)6/6-.-R(x-1)5/5] 

[(x-1  )5/5  +  (Ri-1  )(x-1  )^/4-t-R(x-1  )3/3]exp(-x2/4  R) 

1 

+  |dxx-4[(x  +  1  )5/5-(x+1  )6/6  +  (x-1  )5/5  +  (x-1  )6/6] 

0 

[(x  +  1  )3/3  +  (x-1  )^/4  +  (x- 1)3/3- (x  +  1  )‘^/4]exp(-x2/4  R) 
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oo 

+  Jdxx-4[(x  +  1)5/5-(x  +  1)6/6][(x  +  1)3/3-(x  +  1)4/4]exp(-x2/4R)},  (111) 

1 

where  R=kBT/2m*Ci2,  Ci  is  the  deformation  potential,  and  Nc  is  the  effective  density  of  states 
for  the  conduction  band. 

VIII.  4.  Non-Polar  Optical  Phonon  Scattering 

This  time  one  obtains 

oo 

an  0  ph  =  [87:V2hcooaNch2/3  m  ^''^c^cooK  jdxx^^^ 

■  0 

[(F  +  1  )(x-1)1/20(x-1  )  +  F(x  +  1  )1/2]-4 

[(F+1)2(x-1)(2x-1)0(x-1)  +  F2(x  +  1  )(2x  +  1  )]exp(-hcooX/kBT)} 

oo 

{  jdxx2/2[(F+i  )(X-1  )1/20(X-1  )  +  F(x  +  1  )l/2]-lexp(-hcooX/kBT)}-2,  (112) 

0 

where  F=[exp(hQo/kBT)-1]-i,  and  coq  is  the  optical  phonon  frequency. 

Vlll.  5.  Polar  Optical  Phonon  Scattering 
Proceeding  as  in  Secs.  Vill.2  and  .4,  we  obtain 

oo 

an  0  ph  =  [8jiV2hco|aNch2/3m'5/2c2u|j(  jqxx^ 

[F2(x  +  1  )1''2|n(2x1/2  +  2(x  +  1)1/2) 

+  (F  +  1)2(x-1)l/2|n(2xi/2  +  (x-i)i/2)0(x-i ) ] e x p ( - huix/ksT ) 

[(F  +  1)arcsinh(x-1)i^2  0(x-i)  +  Farcsinh(xi/2)]-4}_  (113) 

Here  coi  is  the  longitudinal  phonon  frequency. 

Vlll.  6.  Intervalley  Scattering 

This  type  of  scattering,  present  in  indirect  bandgap  semiconductors,  transfers  electrons 
from  one  of  the  six  minima  (or  valleys)  of  the  conduction  band  energy  in  k-space  to  one  of  the 
other  five  minima.  Transitions  between  a  valley  and  the  nearest  valley,  which  is  along  the  same 
k-space  direction  in  the  next  copy  of  the  first  Brilloin  zone  in  the  periodic  zone  scheme,  are  of 
the  Umkiapp  type,  and  are  called  g-processes.  Transitions  to  the  four  valleys  present  in  the 
same  zone  along  the  other  two  k-space  directions  are  called  f-processes.  Repeating  a  previous 
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calculation  [10']  on  the  basis  of  the  new  cross-correlation  formula  (106),  we  obtain  for  g- 
processes 

oo 

ttg  =  [8W2hcoijaNch2/3m’5''2c2^^ 

[(F  +  1  )(x-1  )^/20(x-1  )  +  F(x-i-1 

[(F  +  1  )2(x-1)(2x-1  )0{x-1  )-t-F2(x-i-1  )(2x-f-1  )]exp(-hcoijx/kBT)} 

oo 

{  jdxx3/2[(p+i)(x-i)l/20(x-i)  +  F(x  +  1)''^2]-lexp(-hcoijx/kBT)}-2,  (114) 

0 

where  hcojj  is  the  phonon  energy  corresponding  to  the  momentum  difference  required  by  the 
inlervalley  transition.  For  the  corresponding  f-process  we  obtain  [lO'j 

ttf  =  (ko/qo)^«g.  (115) 

where  ko/qo  is  the  ratio  between  the  position  vector  of  a  conduction  band  energy  minimum  in  k 
space,  and  twice  the  distance  of  the  minimum  from  the  Brillouin  zone  boundary,  0.85/0.3  for 
silicon.  There  are  three  g-type  alphas  ogi.  ogz  and  ogs  (from  LA,  TA  and  LO  phonons 
repectively)  and  three  f-type  contributions  afi,af2  and  ao  (from  TA,  LA  and  TO  phonons). 
Their  values  are  given  in  the  graph  below  and  are  a  few  times  larger  than  the  old  values. 

The  various  quantum  1/f  contributions  derived  here  can  be  approximately  superposed  to 
yield  the  resultant  quantum  1/f  coefficient  according  to  the  rule 

OH  =  Ii(|i/pi)2  a.\  (116) 
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IX.  COHERENT  QUAJSTTUM  1/f  CHAOS 

Conventional  quantum  1/f  fluctuations  of  physical  cross  sections  and  process  rates  have 
been  introduced  by  us  as  a  fundamental  infrared  divergence  phenomenon  in  1975  [1].  Some  of 
the  subsequent  publications  [2]-[l3]  have  shown  this  new  effect  to  be  unaffected  by  the 
presence  of  the  thermal  radiation  background  [4],  [5],  some  have  derived  it  with  wave  packets 
[11],  including  a  finite  mean  free  path  [9],  in  second  quantization  [13],  with  the  Keldysh- 
Schwinger  method  [12],  or  in  the  Van  Hove  weak  interaction  limit  [14],  some  derived  its 
characteristic  functional  [7],  or  applied  it  to  the  calculation  [10]  of  mobility  and 
recombination  speed  fluctuations  in  semiconductors  and  semiconductor  devices.  Others  verified 
the  new  effect  experimentally  and  successfully  applied  it  to  electronic  devices. 

The  present  paper  derives  a  related  fundamental  effect  which  we  call  the  coherent  1/f 
effect,  with  elementary  methods  of  quantum  electrodynamics  and  non-relativistic  many-body 
theory.  Our  derivation  uses  the  new  picture  introduced  by  Dollard  [55]  and  generalized  by 
Kulish  and  Fadde'ev  [56]  and  later  by  Zwanziger  [57],  in  agreement  with  earlier  work  by 


Chung  [58]  and  by  Kibble  [59].  in  this  new  picture,  the  asymptotic  Coulomb  interaction  is 
included  in  the  unperturbed  Hamiltonian  rather  than  in  the  perturbation  part.  This  leads  to  a 
more  complex  physical  free  particle  notion  which  includes  a  coherent  photon  cloud,  and  replaces 
the  pole  in  the  propagator  with  a  branch  point.  It  also  leads  to  a  smeared-out  mass  shell.  Using 
this  picture,  we  can  neglect  the  remaining  part  of  the  interaction  if  we  limit  ourselves  to  the 
asymptotic  region  of  large  distances  and  times,  which  are  important  in  the  case  of  1/f  noise. 

For  N  electrons  in  a  Fermi  sphere  shifted  in  momentum  space  by  a  vector  po  and 
occupying  N/2  orbitals  e'P^  the  propagator  derived  by  these  authors  [57]  can  be  reduced  for 
large  time  components  of  x'-x  to  the  non-relativistic  form 

-i«I)o|T\ys'(x')v|-(x)|Oo>  =  5ss'  Gs(x’-x) 

=  (i/V)I{expi[p  (r-r')-p2(t-f)/2m]/fi}np,s 

P 

x{-ip  (r-r‘)/K+i(m2c2  +  p2)l (117) 

Here  a=e2/fic=1/137  is  Sommerfeld's  fine  structure  constant,  np_s  the  number  of  electrons  in 
the  state  of  momentum  p  and  spin  s,  m  the  rest  mass  of  the  fermions,  6ss'  the  Kronecker 
symbol,  c  the  speed  of  light,  x=(r,t)  any  space-time  point  and  V  the  volume  of  a  normalization 
box.  T  is  the  time-ordering  operator  which  orders  the  operators  in  the  order  of  decreasing 
times  from  left  to  right  and  multiplies  the  result  by  (-1)^,  where  P  is  the  parity  of  the 
permutation  required  to  achieve  this  order.  For  equal  times,  T  normal-orders  the  operators, 
i.e.,  for  t=t'  the  left-hand  side  of  Eq.  (117)  is  i<«l>olv|:{x)ys'(x')l'l^o>-  The  state  C>o  of  the  N 
electrons  is  described  by  a  Slater  determinant  of  single-particle  orbitals. 

Consider  first  the  case  t=t'  for  simplicity,  although  only  the  case  of  large  t-t'  can  be 
expected  to  be  experimentally  applicable.  The  pair  correlation  function  can  then  be  decomposed 
as  follows 

<«1’oivl(>:)vl'’(x’)vs’(x')Vs{x)l‘I>o>  =  <‘1>olv^(x)vs(x)|Oo><'l>olv^'(x')Vs'(x')l^o> 

-  <‘f’olvl(x)Vs'(x')|C)o><Oolv^'(x')v)/s(x)|cI>o>.  (118) 

The  first  term  can  be  expressed  in  terms  of  the  particle  density  of  spin  s,  n/2  =  N/2V  = 
<'f’olv^(x)vs(x)|‘&o>.  while  the  second  term  can  be  expressed  in  terms  of  the  Green  function 
(117)  in  the  form 

Ass'(x-x')  =«:>olv^(x)vJ-(x')vs'(x')Vs(x)|Oo>  =(n/2)2  +5ss'  G s(x'-x)Gs(x-x') .  (119) 
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The  "relative"  autocorrelation  function  A(x-x’)  describing  the  normalized  pair  correlation 
independent  of  spin  is  obtained  by  dividing  by  n^  and  summing  over  s  and  s' 

A(x-x’)  =  1  +  (1/n2)XGs(x-x')Gs(x'-x) 

s 

=  1  -  (1/N2)X  Z{exp[i(p-p')(r-r')/fi]}np,snp*.s 

s  pp' 

X  Ip  (r-r')/Kl“/’'lp’(r-r')/K|«/’'.  (120) 

Here  we  have  used  Eq.  (117).  The  low-wavenumber  part  A\  of  this  relative  density 
autocorrelation  function  is  given  by  the  terms  with  p=p'. 

Ai(x-x’)  =  1  -  (1/n2)X  Inp,slp(r-r')/Kl2<x/n  (121) 

s  p 

=  1  -  (2/N2)[V/(27ih')3]  |d3p|(p  +  po)(r-r')/Ii|2a/7t 

P<Pf 

a  a  a 

=  1  -  (2/N2)[V/(2rttt)3]  Jdpi  |dp2  Jcip3l(P3  +  Po3)(r-r')/h'|2“/" 

-a  -a  -a 

=  1  -  (2/N2)[V/(27th')3]8a3|a(r-r')/h'|2“/rt 

=  1  -  (1/N)[(::/6)‘'^PF|r-r'l/K]2“/’t.  for  pF»po3: 

and 

=  1  -  (1/N)|po(r-r')/K|2“'»'  forpF«po3: 

In  view  of  the  smallness  of  2a  n,  for  Pf>>Po3  we  have  integrated  in  cartesian  coordinates, 
approximating  the  Fermi  sphere  by  a  cube  of  side  2a  with  a=  (7r/6)'''2pF.  The  result  is 
practically  independent  of  Po  and  of  Po3=Po('’-f’)^l*’-''’l-  '•’he  factor  (7i/6)2“''3rt  can  be 
neglected.  For  Pf«Po3  we  used  the  mean  value  theorem  for  estimating  the  integral  over  d^p  in 
spherical  coordinates.  In  all  cases  the  autocorrelation  decreases  very  slowly  from  1  when  |r- 
r'l  is  increased  to  very  large  values.  Writing  the  rectangular  bracket  in  Eq.  (122-123)  as  an 
exponential  function  of  its  logarithm,  expanding  the  resulting  exponential,  and  keeping  only  the 
first  term,  we  obtain  with  pF/HskF 

Ai(x-x')  =  (1-2/N)  +  (6/;c)2«'3''{[1/kF|r-r'l]2«/:t}/N 

oo 

=  (1-2/N)  +  (6/7t)2a'37t(i/kp]2a/7t(2a/N;i)|cos[k|r-r'|]dk/kl-2a/7t. 


(122) 

(123) 
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oo 


={N-2  +(2a/7:)J[k/kF]2«/"cos[k|r-r'|]dk/k}/N 


for  Pf>>Po3I 


(124; 


={N-2  +(2a/7i:) J[k/ko]^“''^cos[kpo(r-r’)/po]  dk/kj/N  for  Pf<<Po3:  (125) 

0 

Here  we  have  used  a  well-known  Fourier  integral  [60]  and  we  have  introduced  ko^poa/K. 
According  to  the  Wiener-Khintchine  theorem,  the  coefficient  of  the  cos  gives  the  spectral 
density.  To  get  it  for  the  fractional  fluctuations  Sn/n,  we  divide  by  the  constant  term  N-2 

S5n/n(k)  =  [2a/7ik(N-2)][k/K]2«/:t^  (126) 

where  K=kF  for  Pf»Po3  as  in  the  case  of  metals  with  spherical  wave  symmetry,  and  K=ko  for 
Pf<<Po3-  Although  inapplicable,  this  pure  1/k  spectrum  is  the  wave-number  equivalent  of  the 
coherent  quantum  1/f  noise  derived  earlier  [61],  [62]  in  excellent  agreement  with  the 
experiments  on  large  electronic  devices  [46],  Due  to  2aA«1  the  second  factor  is  practically 
unity  and  of  no  importance,  except  for  eliminating  the  logarithmic  divergence  from  the  spectral 
integral.  This  wave  number  spectrum  also  entails  a  1/f  frequency  spectrum  obtained  by 
writing  dk/k=df/f,  as  was  shown  in  detail  in  a  previous  paper  [13].  For  equal  times  our  result 
can  not  be  expected  to  be  valid,  due  to  the  asymptotic  character  of  Eq.  (117).  We  shall  now 
derive  the  1/f  spectrum  directly  below. 

If  t?it',  Eq.  (118)  is  replaced  by 

<4>oITvI(x)Vs(x)v^'(x')vs'(x')14>o>  =  «t’olvl'(x)Vs(x)|<l>o><4>ok|’(x')\|/s'(x')|4)o> 

-  <Oo|Tvs'(x')v^(x)l<^o><‘l>olTvs(x)y|-(x’)l(I)o>.  (127) 

Eq.  (119)  remains  the  same,  except  for  the  middle  part  which  is  replaced  by  the  left  hand  side 
of  Eq.  (127).  Eq.  (120)  becomes  now 

A(x-x’)  =  1  -  (1/n2)V  Gs(x-x')Gs{x'-x) 

s 

=  1  -  (1/n2)X  I{expi[(p-p')(r-r’)-(p2-p'2)(t-f)/2m]/Fi}np,snp',s 

s  pp' 

x{p(r-r')/h'-(m2c2+p2)1  /2(j.i-)(c/h')}a/7i 

x{p  ■(r-r’)/h'-(m2c2  +  p'2)1/2(t.f)(c/h')}a/Jt.  (128) 
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Here  we  have  used  again  Eq.  (117).  The  low-frequency  and  low-wavenumber  part  A|  of  this 
relative  density  autocorrelation  function  is  also  given  by  the  terms  with  p=p'. 

Ai(x-x')  =  1  -  (1/n2)X  Snp,s 

s  p 

x{p  (r-r')/h'-(m2c2  +  p2)1/2(|.t')(c/jf)j2a/7i  (129) 

=  1  -(2/N2)[V/(27:K)3  Jd3p|(p  +  po)(r-r‘)/fi-[m2c2-H(p  +  po)2]i/2(t-t')(c/K)i2a/7t 

fXPF 

=  1  -  (1/N)lp  o(r-r')/K-mc2T/h'pa/it  for  pp«|po3-mc2T/z|.  (130) 

Here  we  have  used  the  mean  value  theorem,  considering  the  laJ-n  power  as  a  slowly  varying 
function  of  p  and  neglecting  po  in  the  coefficient  of  t  =t-t’,  with  z=|r-r'|.  Writing  the  power 
again  as  an  exponential  function  of  its  logarithm,  expanding  the  resulting  exponential,  and 
keeping  only  the  first  term,  we  obtain  with  (h7mc2)|po(r-r’)/K-mc2x/h'|  =  6  =  |x-po(r- 
r')/mc2| 

Ai(x-x')  =  1-[(mc2/h')e]2“/’t/N 

oo 

=  (1-2yN)  +  [K/mc2]2a/n(2a/N7t)|coS[a)6]dco/col“2a/7t_ 

0 

oo 

=■  {N-2  +  (2a/7t)J'[h'c}/mc2]2«/JtcQ5|(j30j  dto/ajj/N.  (131) 


According  to  the  Wiener-Khintchine  theorem,  the  coefficient  of  the  cos  gives  the  spectral 
density.  To  get  it  for  the  fractional  fluctuations  5n/n,  we  divide  by  the  constant  term  N-2 

S5n/n(k)  =  [2a/7ia)(N-2)l[h'(o/mc2]2a/7t  (132) 

This  result  is  also  applicable  for  the  particle  current  fluctuation  spectrum.  Indeed,  for  current 
density  fluctuations  5j  we  include  a  (fi/mi)V  in  front  of  each  of  the  two  y  operators  in  Eq. 

(127),  a  factor  ppVpo^  in  Eq.  (128)  after  the  summation  signs,  a  factor  (p/po)^  in  the  first 
form  of  Eq.  (129),  a  factor  (p+Po)^/Po^  in  the  second  form,  and  no  changes  in  Eqs.  (130)- 
(131).  Eq.  (132)  becomes 

[2tt/7m(fii2j.IIfii^mc2pa/»i.  (13  3) 
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This  result  coincides  with  our  earlier  theoretical  result  for  coherent  quantum  1/f  noise  if  we 
replace  N  with  N-2.  The  validity  of  this  equation  is  restricted  to  low  frequencies  and  wave- 
numbers.  This  equation  is  in  excellent  agreement  with  mobility  and  diffusion  1/f  noise  in  large 
devices. 

Finally,  we  consider  the  errors  caused  by  the  neglect  of  higher  order  terms  in  the 
expansion  of  the  exponential  functions  resulting  from  Eqs.  (122)-(123)  and  (130).  For  a 
thermal  eleotron  and  r=1cm  in  Eq.  (6')  we  get  kr=10''2  and  (2a/n)ln(kr)=0.12,  yielding  an 
error  of  12%.  For  t=10''^s  in  Eq.  (130),  which  is  the  age  of  the  universe,  we  get  an  error  of 
40%  as  the  upper  limit.  Therefore  a  more  exact  treatment  is  of  some  interest.  Using  the 
identity  [60] 


oo 

e2a/n  ^  [-(2a/7c)  |a)-2“/^COS(0co)d(o/(D] 

COo 

x{cosa  +  (2a/;t)^(e(0o)2"-2«/"[(2n)!(2n-2a/7r)]-1}■^  (134) 

n=0 

with  arbitrarily  small  cutoff  coq.  we  obtain  from  Eq.  (130)  the  exact  form 


A|(x-x')  =  1  +  [(2a/;cN)  J(mc2/K(o)2“/"cos(eco)dco/a)] 

“o 

OQ 

x{cosa  +  (2a/K)2^(ea)o)2"-2«'"[(2n)l(2n-2a/K)]-1}-^  (135) 

n  =  0 


This  would  indicate  a  ci)''''2'^''  spectrum  and  a  1/N  dependence  of  the  spectrum  of  fractional  n 
and  j  fluctuations,  if  we  neglect  the  curly  bracket  in  the  denominator  which  is  close  to  unity  for 
very  small  wq-  We  thus  realize  that  the  unusual  N-2  dependence  in  Eqs.  (132)-(133)  is 
caused  by  the  forced  introduction  of  the  integrable  spectrum  in  place  of  the  a)"'“2a/7i 

spectrum.  Due  to  the  smallness  of  a  both  forms  coincide  in  practical  applications.  Eq.  (133) 
for  the  coherent  QED  chaos  process  in  electric  currents  can  thus  be  written  also  in  the  form 

Ssi/i(k)  ==  [2a/7ta)N1[mc2/R0j]^n'>t  -  2a/7ia)N  =  0.00465/o)N.  (13  6) 

This  result  derived  directly  earlier  [61],  [62],  is  in  excellent  agreement  with  the 
measurements  [46],  in  large  [see  [62]  for  a  definition  of  "large"  or  "extended",  and  for  an 
interpolation  with  the  conventional  quantum  1/f  chaos  effect}  devices  such  as  large  n+p  Hg^. 
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xCdxTe  infrared  detector  diodes.  It  is  also  close  to  the  empirical  value  of  0.002/coN  observed 
earlier  by  Hooge  [63]  in  semiconduotors  and  metals.  Being  observed  in  the  presence  of  a 
constant  applied  field,  these  fundamental  quantum  current  fluctuations  are  usually  interpreted 
as  mobility  fluctuations. 
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